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Abstract 



(N 
> 

00 ■ We derive the leading asymptotic behavior and build a new series representation for the Fredholm determinant 

in I of integrable integral operators appearing in the representation of the time and distance dependent correlation 

functions of integrable models described by a six-vertex /?-matrix. This series representation opens a systematic 
way for the computation of the long-time, long-distance asymptotic expansion for the correlation functions of the 
aforementioned integrable models away from their free fermion point. Our method builds on a Riemann-Hilbert 
based analysis. 



1 Introduction 

Highly structured determinants appear as a natural way for representing the correlation functions in integrable 
models that are equivalent to the so-called free fermions. It was already shown by Kaufman and Onsager that cer- 
tain two-point functions of the 2D-Ising model can be represented by Toeplitz determinants [30]. Then Montroll, 
Potts and Ward [44] made this observation more systematic by expressing the so-called row-to-row two-point 
function of this model in terms of a Toeplitz determinant. It was observed by Lieb, Mattis, Schultz [40] that such 
Toeplitz determinant-based representations also hold for the so-called XY model. Then, the systematic study of 
the correlation functions of the impenetrable Bose gas, the XY model or its isotrpoic version the XX model lead 
to the representation of various correlators in terms of Fredholm determinants (or their minors) of the so-called 
integrable operators [7, 8, 36, 39, 42, 47]. Such types of Fredholm determinants also appear in other branches of 
mathematical physics. For instance, the determinant of the so-called sine-kernel acting on an interval J is directly 
related to the gap probability (probability that in the bulk scaling limit a given matrix has no eigenvalues lying in 
/) in the Gaussian unitary ensemble [21]. Integrable integral operators [10] are operators of the type I + V where 
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the integral kernel V takes a very specific form. This fact allows for a relatively simple characterization of the 
resolvent kernel and often for a construction of a system of partial differential equations satisfied by the associated 
Fredholm determinant or minors thereof [14, 24, 26, 29, 49]. 

In all of the aforementioned examples, the integrable integral operators / + y act on some curve with a 
kernel V {A, fi) depending, in an oscillatory way, on a parameter x. In the previous examples a lot of interesting 
information can be drawn out of the asymptotic behavior of det [/ + V] for large values of x. For instance, when 
deaUng with the correlation functions of integrable models, x plays the role of a spacial and/or temporal separation 
between the two operators entering in the correlation function. In such a case, computing the large-x asymptotic 
expansion of the associated Fredholm determinants, allows one to test the predictions of conformal field theories. 
The form of the asymptotic behavior of the pure sine kernel determinant log det [/ + 5 ] was strongly argued in 
[6, 19] and then proven, to some extend, using operator methods [3, 20, 51]. Also, the discovery of non-Unear 
differential equations of Painleve V type for this determinant [29] allowed to compute many terms in the large-x 
asymptotic expansion of the associated correlation functions [29, 41, 42, 43] 

However, a really systematic and efficient approach to the asymptotic analysis of various quantities related 
to integrable integral operators I + V has been made possible thanks to the results obtained in [26]. There it 
was shown that the analysis of such operators can be reduced to a resolution of an associated Riemann-Hilbert 
problem (RHP). The jump contour in this RHP coincides with the one on which the integral operator acts and the 
jump matrix is built out of the functions entering in the description of the kernel. In this way, one deals with a 
RHP depending on x in an oscillatory way. The asymptotic analysis of their solutions is possible thanks to the 
non-linear steepest descent method of Deift-Zhou [16, 17]. It is in this context that the full characterization of the 
leading asymptotic behavior of Fredholm determinants of kernels related to correlation functions in free-fermion 
equivalent models (the long-distance, long-time/long-distance at zero and also non-zero temperature) has been 
carried out in the series of papers [4, 9, 23, 25, 27, 28]. Also, it was shown that the value of Dyson's constant 
arizing in the asymptotics of the pure sine kernel determinant can be obtained in the RHP-based framework as 
well [13, 38]. 

This paper deals with an extension of these analysis to the case of a Fredholm determinant of an integrable in- 
tegral operator whose integral kernel has a more involved structure then in the aforementioned cases. We call our 
kernel the time-dependent generalized sine kernel. The Fredholm determinant we analyze arises in the represen- 
tation of the zero temperature long-distance/long-time asymptotic behavior of two-point functions in a wide class 
of integrable models away from their free fermion point. In particular, its asymptotics expansion (and especially 
the new series representation that we obtain for it) plays a crucial role in the computation of the long-time/long- 
distance asymptotic behavior of these two-point functions. Therefore the results in this paper can be seen as a first 
step towards such an asymptotic analysis of two-point functions. Let us be slightly more specific. 

In a wide class of algebraic Bethe Ansatz solvable models, one is able to compute the so-called form factors 
(matrix elements of local operators) and represent them as finite-size determinants [35, 46] . It has been shown in 
[36] that, for free fermion equivalent models, it is possible to build on these representations so as to explicitly sum- 
up the form factor expansion and compute the zero-temperature (and even the non-zero temperature) correlation 
functions of the model. In the limit of infinite lattice sizes, a two-point function is then represented by a Fredholm 
determinant (or its minors) of an integrable integral operator I + V acting on some contour 'rf determined by the 
properties of the model. For time and space translation invariant models, the kernel V depends on the distance 
separating the two operators as well as on the difference of time evolution between them. One can show that for 
general free-fermion type models, the integral operator I + V associated with the form factor expansion of the time 
and space dependent two-point functions acts on a finite subinterval [ ; ^ ] of R and its kernel V belongs to the 
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class of kernels 

sin Ittv (/()] sin I nv (li) | 

(1.1) V(A,fi)^4 \] ' \ '^'^ {E{A)e{M)-E(ju)e{A)} . 

2m {A — n) 

There v is some function encoding the fine structure of the excitations above the ground state whereas e and E are 
oscillating factors. The function E is expressed in terms of e 

(1.2) E{A) = ie{A){ I ^tlS!l + tlS^coi[nv{A)V, 



'1/ 



2n s-A 



The functions v and e just as the integration curve '^e appearing in (1.2) depend on the specific model that one 
considers. We will give more precision about their properties in the core of the paper. We stress that, for free- 
fermion equivalent models, v{A) is some constant and e takes a simple form. It was in such a context that the 
asymptotic analysis of det [/ + V] has been carried out previously. 

When considering integrable models that are away of their free-fermion point, as it has been shown in [34, 37], 
it is as well possible to build on the finite-size determinant representation for the form factors of local operators in 
integrable models out of their free fermion point so as to sum up the form factor series over the relevant sector of 
excited states. Even though Fredholm determinant based representations are no longer possible in models away 
from their free fermion point, the above procedure leads to a representation for the time and space dependent 
two-point functions in terms of series of multiple integrals which take the generic form 



n>0 „ a=l Z, a=l ^ a=l a=\ 



1 
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There, and are some contours in C surrounding [ ; w is a holomorphic function in a neighborhood 
of [ ; ^ ] U and fn some function. In fact, an analogous type of series representation (1.3) have been first 
obtained, through other method, in [33], this in the case of the space dependent functions only. That paper also 
proposed an approach to the extraction of the large- asymptotic behavior out of the "pure-distance" analogue of 
the series (1.3). 

An important observation made in [34, 37] is that the intermediate computations can be shown to boil down to 
effective generahzed free fermionic models. As such, they involve, again, the Fredholm determinants of operators 
/ -I- y with V given by (1.1). However, then, the functions v and e become much more complex that at the free 
fermion point. In some sense, the approach of [34, 37] shows that kernels (1.1) appear as a natural basis of 
special functions allowing one to represent the correlation functions of a wide class of interacting {ie away from 
their free fermion point) integrable models as certain (infinite) Unear combinations thereof. Therefore, the main 
motivation for our study of the time-dependent generalized sine kernel (1.1) is to obtain a convenient and effective 
representation - that we call the Natte series^ - for the associated Fredholm determinant. The Natte series allows 
one to re-sum the aforementioned decompsition of the representation (1.3) as a linear combination of Fredholm 
determinants into some compact and explicit form. The latter provides one with a new type of series of multiple 
integrals representation for two-point functions. Moreover, it is built in such a way that one is able to read-off the 
asymptotic behavior of the correlators out of it. Therefore, the results established in this paper provide one of the 
fundamental tools that are necessary for carrying out the long-distance and large-time asymptotic analysis of the 
two-point functions in massless integrable models proposed in [34, 37]. 



'The origin of this name issues from the so-called pig-tail (or braid) hairstyle that is called Natte in French. A braid is a specifically 
ordered reorganization of the loose hair-do style. Similarly, the Natte series reorganizes the Fredholm series in a very specific way, so that 
the resulting representation is perfectly fit for carrying out an asymptotic expansion. 
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This article contains two main results. We first derive the leading asymptotic behavior of the Fredhohn de- 
terminant of / + V understood as acting on {[—q;q]), with q < +00. This sets the ground for the second 
main result of the article. Namely, we derive a new series representation -the Natte series- for the Fredholm 
determinant. This series is converging rather fast in the asymptotic x +00 regime. Its main advantage is to 
provide a direct (ie without the need to perform any additional analysis) approach to the asymptotic expansion of 
the determinant. As already stressed out, this series representation plays a crucial role in the computation of the 
lagre-time/long-distance asymptotic expansion of the two-point functions in integrable models corresponding to 
a six-vertex /?-matrix. Also, the very form of the asymptotic expansion stemming from the Natte series proves 
several conjectures relative to the structure of the asymptotic expansions for certain particularizations (for spe- 
cific values of y, and e) of such Fredholm determinants [45, 48]. Also, upon specialization, it yields the general 
structure of the asymptotic expansion of the fifth Painleve transcendent associated to the pure sine kernel [14, 29] 

This article is organized as follows. In section 2, we outline the main assumptions that we rely upon throughout 
the article and give a discussion of the class of functions e that we deal with. After introducing several notations, 
we present the two main results of the paper. The remaining part is of technical nature. In section 3, we present 
the RHP problem that is at the base of the asymptotic analysis of det [/ -1- V] and the construction of its Natte 
series. We also outhne the chain of transformations corresponding to the implementation of the Deift-Zhou [17] 
steepest-descent method. In section 5, we build the various local parametrices. This brings the original RHP 
into one that can be solved through a series expansion of the associated singular integral equation [5]. The latter 
naturally provides the large-x asymptotic expansion of the solution. We build on these results so as to derive 
the leading asymptotic expansion of the Fredholm determinant in section 6. Finally, section 7 is devoted to the 
construction of Natte series for the Fredholm determinant of I-i-V. In particular, we establish the main properties 
of such series. We then give a conclusion and discuss the further possible applications. In appendix A, we recall 
all the properties of the special functions that we use in this article. In appendix B, we gather some proofs relative 
to the structure of the large-x asymptotic expansions of certain matrix valued Neumann series representing the 
solution to a singular integral equation of interest to us. In appendix C, we establish bounds for certain matrices 
appearing in our analysis. 



2 The main results and assumptions 



In this article, we will focus on the case where the function e takes the form 



(2.1) e-'{A) = e'— 



:MA) . g(A} 
2 



+ - 



e (A) is quickly oscillating in the x — > -i-cxd limit and the function g entering in the definition of e (A.) has been 
introduced so as to allow for some finite, x-independent oscillatory behavior of the function e (A). The principal 
value integral apperaing in the definition of £ (1.2) is carried out along a curve ^£ which corresponds to a sUght 
deformation of the real axis and is depicted in Fig. 1. Under the forthcoming hypothesis, such a contour allows to 
strengthen the convergence of the integral defining E at infinity (in the case of K., the convergence would be the 
one of an oscillating non-absolutely integrable power-law whereas it is exponentially fast along "^g). 



2.1 The main assumptions 



Throughout this paper, we make several assumptions on the function u, g and v entering in the description of the 
integrable kernel (I.I). 

• There exists an open neighborhood UofR such that u and g are simultaneously holomorphic on U. 
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• The function g is bounded on U. 

• The function u is real valued on R and has a unique saddle-point in U located at Aq e R. This saddle-point 
is a zero of u' with multiplicity one, that is to say 3 ! /Iq e R : u'(Ao) - and u"(Ao) < 0. We also assume 
that the saddle-point lies away from the boundaries: Aq + ±q. 

• M is such that, given any t] > 0, e"'"^'^^ decays exponentially fast in A when ±3 (A) > 6 > for any fixed 

5 > 0, and 21 (/I) ^ +oo. 

AeU 

• The function v is holomorphic on U and such that sin [nv (A)] has no zeroes in some open neighborhood of 

[-q;q] lying in U. 

• The function v has a "sufficiently" small real part at ±q, ie \% [v (±^)]| < 1 /2. 

For technical reasons, one has to distinguish between two situations when the saddle point Aq is inside of 
] -(7 ; <j [ or outside. Following the tradition we refer to the first case {-q < Ao < q) as the time-like regime and to 
the second one {\Aq\ > q) as the space-like regime. Actually, in this article we will only consider the case where 
Aq> q. Also, we do not treat the Umiting case when Aq - ±qd& this would require a significant modification of 
our approach. 



2.2 The main result 



We now gather the main results of this paper into two theorems. 

Theorem 2.1. Let V{A,ii) be as in (1.1) and I + V act on i\^—q;q\}. Then, under the assumption stated in 
section 2.1, the leading x — > +oo asymptotic behavior o/det [/ + V] reads: 



(2.2) det [/ -H y] e 
{-<l^'<l^ 



^BAv,u\\\+0 



logx\\ bi [v,u,g] 



3 

X2 



B, [y,M] l + O 



logx] 



^ ^i4ui<!)-u(-ci)]+g(q)-g(-g)g^ + 1, m] |l + O |i^JJ + e'4"(-^)-"(^)]+^(-9)-«Ws^ [y - 1, «] |l + O {^■^''^ 

The functional [v, u] takes the form 

(2.3) [V, u] = e^'M GHi + v (g)) GHi - v j-g)) .Hi^yA-,)) 

{2gx{u'{g) + jO+)]''^^^ {2gxu'{-g)Y^-'i^ 



It is expressed in terms of the Barnes G function [1] and the auxiliary functional 



(2.4) Ci [v] = - I d^dju 

-1 



<4 



v'{A)v{jD-v' (pi)v{A) 
A- jj. 



+ v{g) 



-v{A) 



&A + V i-g) 



Cvi-qy- 
J q + 



-q)-v{A) 



A 



di 



-1 



The functional b\ [y, u, g\ takes different forms depending whether one is in the so-called space-like regime (Aq > g) 
or in the time-like regime (A() e]-g■,g[): 



(2.5) bi [v, u, g\ 



1 



v{-g) So 



vig) 



^J-2nu"{Ao) 



u'(-q)(Ao+qfS^ u'ig)iAo-gfS+ 
v{-q) <$- _ v{q) 
, u'{-q){AQ + qf Sq u' (q) (Aq - qf So 



time - fike 
space - Uke 
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There, we agree upon 



(2.6) . [2,x«'(,) . .-O^]-^^) e\,) (e--^^) - l) ^^i^ exp | - 2 f'^ ^^Mzl^d,} , 

^ 'T{l+v{q)) \ J^q q-n } 

ni, ^ (e-^'>"-(-^) - 1) ra+v(^ r v(-^)-v(/^) , 1 

(2.7) 5- = T-;^e (-f?) ITT. : — :t exp < 2 du \ , 



(2.8) 5o = . (^)e.(^-— -^j exp|-2j_^— ^— — d^xJ 



space - like 



The proof of this theorem will be given in section 6.2. It heavily relies on the asymptotic analysis of the RHP 
associated with V that will be carried out in sections 3 to 5. 

Above, the iO^ regularization is important only in the time-like regime as then u'{q) < 0. It allows one for 
a non-ambiguous definition of the power-laws appering above. In the space-Uke regime, the iO^ regularization 
makes no difference . 

A special limit of the kernel (1.1) can be related to the generalized sine kernel studied in [32]. Indeed, when 
the saddle point is send to infinity, by deforming slightly the contours 'Ioe, the function E can be seen to 
be proportional to e ^ up to corrections that are uniformly 0{x'°°) on [-^;^]. In particular, one has that the 
X — > +00 asymptotic expansion of the two Fredholm determinants coincide in this limit. This can be seen directly 
by inspection of our formulae, at least in what concerns the leading asymptotics. 

A specific case of our kernel u{X) = A- tA? jx, g = 0, q - \ and v = est has been studied in the Utterature in the 
context of its relation with the impenetrable fermion gas [4]. Upon such a specialization, our results agree with 
the coefficients of the asymptotic expansion obtained in that paper. 

The second main result obtained in this paper is the Natte series representation for the Fredholm determinant. 

Theorem 2.2. Under the assumptions stated in section 2.1, the Fredholm determinant of the operator I + V where 
the kernel V is given by (l.l) admits the below absolutely convergent Natte series expansion. In other words, there 
exists functionals fin [v, e^, e""j such that 

(2.9) det [/ + V] [y, u, g\ = det [/ + V]^^^ [v, |l + 2 t^' "j| • 
There 

(2. 10) det [/ + Vi^^ [v, u, g] = [v, m] • exp | [ixu'(A) + g'(A)] v (A) d^j . 

A more detailed structure of the functionals 'Hn can be found in the core of the text, formulae (7.10). One has the 
following estimates for the functionals YHn [v, e^, .ic]| < [m {x)Y, with m{x) = (x^*) being n-independent and 

(2.11) w = ^min(l/2,I -w-2max|!Ry(e^)|) with w = 2sup{|!R [y(i) - y(e^)]| : \A- eq\^ 5 , ±] , 



where 6 > is taken small enough. Hence, the series is convergent for x large enough. 
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The junctionals [v, e^, u] take the form 

L2J ixm\u(q)—u(—q)\ 

(2.12) 7/„[v,e^«]=^r[v,e^«]. 2 o4v(.Hv(-.)] ^"'"^ 

'"=-[!] 

ZjZj Zj ^2m[v(?)+v(-g)] x2r,(b-p)y(-q)-2npy(q) ^" Lnc,MJ. 

^,= 1 p=0m=;,-[5] ' 

Above, we agree upon r} = \ in the space-like regime and rj = -\ in the time-like regime. There fi^^ [v, e*, x] — 
O {x~°°) and the junctionals 'H^^ [v, e^, u] and "T/^™'^'*-* [y, e^, u] admit the asymptote expansions 

[v, u] ~ Y -K^r; [y, u] with <W<r; [y, «] = O i^Hi^^; , 

//I \n+r-2(m+^)) \ 

(2.13) <H^'*'^) [y, e^, «] ~ ^ ^^r/''"^ [^'^ «] with -H^r/'^^) [y, u]=0\ ^ ^ . 



This theorem, together with a more expUcit expressions for the functionals TY^, will be proven in section 7. 
Here, we would however like to comment on the form of the asymptotic expansion. Indeed the above asymptotic 
expansion is not of the type usually encountered for higher transcendental functions. In fact, the large x-behavior 
of the functionals "H,, (y, e^, u) and hence of the determinant det [/ + V] contains a tower of different fractional 
powers of x, each appearing with its own oscillating pre-factor. Once that one has fixed a given phase factor and its 
associated fractional power of x, then the corresponding functional coefficients [y, e^, w] or <j-(j^'^'P^ [y, e^, u] 
admit an asymptotic expansion in the more-or-less standard sense. That is to say, each of their entries admits an 
asymptotic expansion into a series whose r* term can be written as Pr+n (log x) /x"'*''' with P^+„ being a polynomial 
of degree at most r + n. One of the consequences of such a structure is that an oscillating term that appears in a 
sense "farther" (large values of n) in the asymptotic series might be dominant in respect to a non-oscillating term 
present in the "lower" orders of the Natte series. This structure of the asymptotic expansion proves the conjectures 
raised in [45, 48] for certain specializations of this kenrel. Also, upon specialization, it yields the general structure 
of the asymptotic expansion of the fifth Painleve transcendent associated to the pure sine kernel [29]. 

The series representation (2.9) might appear abstract since there is no generic simple expression for the func- 
tionals "H„, n > 1. However, the slightly more explicit (but also more cumbersome so that we did not present it a 
this point) characterization of the functionals 'Hn, gives a thorough and explicit description of the way "74 acts on 
e^. This characterization, together with the overall form of the Natte series (2.9), is enough to build a multidimen- 
sional deformation of (2.9) which describes a class of correlation functions appearing in integrable models away 
from their free fermion points. The very fact that the series representation one starts with has good properties 
from the point of view of an asymptotic analysis (for instance it immediately provides the leading asymptotics) 
leads to a multidimensional deformation which has basically the same good properties in respect to the asymptotic 
analysis, in the sense that it admits an expansion of the type (2.9), (2.12), (2.13). As a consequence, the long- 
time/long-distance asymptotic behavior of two-point functions in an interacting integrable model can be simply 
read-off by looking at the multidimensional series. 



2.3 Notations 



We now introduce several notations that we use throughout the article. 
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• £)zo,(5 = {z e C : |z - zol < <5} is the open disk of radius 6 centered at zq. dD^^^s stands for its canonically 
oriented boundary and -dD^^s for the boundary equiped with the opposite orientation. 

• 0-3, cr* and I2 stand for the below matrices 

'-^> ^^-(Vo) ^--{ro) '-(j^)- 

• Given an oriented curve in C, F C^) stands for a small counterclockwise loop around the curve This 
loop is always chosen in such a way that the only potential singularities of the integrand inside of the loop 
are located on ^. For instance, if consists of one point A, then F (^) can be taken as dDx,6^ for some 
(5 > and small enough. 

• When no confusion is possible, the variable dependence will be omitted, ie u {A) = u, g (A) = g, etc. 

• log refers to the ] -:/r ; ;r [ determination of the logarithm, and it is this determination that is used for defining 
powers. 

o 

• Given a set U, U refers to its interior and U to its closure. 

• H+, resp. EI_, stands for the upper |z e C : 3 (z) > 0|, resp. lower |z e C : 3 (z) < 0}, half -planes. 

• Given matrix valued functions M{A), N{A), the relation M{A) - O {N(A)) is to be understood entry-wise 

MkfiA) = OiNkeiA)). 

• Given an oriented curve "rf, one defines its -1- (resp. -) side as the one lying to the left (resp. right) when 
moving along the curve. Above and in the following, given any function or matrix function /, /+ (A) stands 
for the non-tangential hmit of / (z) when z approches the point Ae'i^ from the ± side of the oriented curve 

• Given a piecewise smooth curve ^ and matrix M with entries in Crf), p = 1, 2, cx5, we use the canonical 
matrix norms (t stands for Hermitian conjugation): 

(2.15) ||M||i~(^) = niax ||Miy||^„^^^ , l|M||i2(^) = J|tr [MtM]||^i^^^ and ||M||ii(^) = niax ||m,-^H^i^^^ . 

• The distance between any two subsets A, S of C will be denoted by d(A, B) = mf[\x-y\ : x e A , y e B}. 

2.4 Several remarks 

It now seems to be a good place so as to gather several remarks in respect to our assumptions. 

• The assumptions on the type of the saddle-point at Aq guarantee that there exists a local parametrization for 
u(A) around Aq, u(A) - u(Ao) - -a/(A) with cj(A) - (A- Ao)h (A), where h (Aq) ^ and h is holomorphic on 
0^0,6 for some 6 > 0. 

• As it will become apparent from our asymptotic analysis, given functions u, v, g satisfying to all the hypoth- 
esis, one has that det [I + V] for x large enough. 

• The assumption on the number of saddle-points and their order can be relaxed in principle. RHP with 
multiple saddle-points have been considered in [50]. This work was later extended to the case of less regular 
functions and higer order saddle-points in [18]. 
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• The restriction on the real part of v in the vicinity of ±q is of technical nature. It allows us to avoid the 
analysis related to the so-called ambiguous Fisher-Hartwig symbols. The method for dealing with such 
kinds of problems in the framework of Riemann-Hilbert problems has been proposed in [11, 12]. The cases 
where 'K (v (±^)) > 1 /2 could in principle be treated along these techniques, but we chose not to venture 
into these technicalities. 

• We have depicted the contour "^e appearing in principal value integral in (1.2) on Fig. 1. This contour is 

chosen in such a way that the integral is converging exponentially fast at infinity. This avoids us unnecessary 
complications and corresponds to most, if not all, situations that can arize in interacting integrable models. 

• In the case of kernels involved in the representation of the two-point functions in integrable models, the 
function u takes the form u{A) - p (A) - ts (/I) /x. p corresponds to the momentum of excitations whereas e 
corresponds to their energy. The parameter t plays the role of the time-shift between the two operators and x 
that of their distance of separation. In general, one is interested in the large-distance/long-time behavior of 
the two-point function in the case where the ratio tjx is fixed. In such a limit, for many models of interest, 
the function u has a unique saddle-point on R. This physical interpretation can be seen as a motivation for 
certain of our assumptions. 

• It is not a problem to carry out the same analysis in the case where the contour given in Fig.l is replaced 
by 'W''^^ ^ n {z € C : |2l (z)| < w}, with w € R+ such that q, -q and belong to ] -w ; w [. Up to minor 
modifications due to such a truncation of the remote part of the contour, the results remain unchanged. 



3 The initial RHP and some transformation 
3.1 The RHP for 

The kernel of any integrable integral operator can be recast in a form allowing one to give a convenient character- 
ization of the kernel R {A, n) of the resolvent operator I -Riol + V. 

Namely, in the case of the kernel V given in (1.1) one sets 




Figure 1. Contour '^e for the definition of E. 



(3.1) 




so that the kernel V is expressed as the scalar product: 



(3.2) 



V{A,n)^ 



{EHA)\E^(j,)} 
A- 
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The resolvent I — Rof I —V exists if det [/ + V] ^ 0. In that case, one defines | (A) ) as the unique solution to 
the integral equation: 

(3.3) iF'^iM))^ I V(M,MF\A))dA^\E-(M)) , | ) , ^ sin [.v (i)] U) j 

where the integration is to be understood entry-wise. ^ (A) | corresponds to the solution of the integral equation 

where | F^ (A) ^ has been replaced with ^ (A) |. It was shown in [26] that the resolvent kernel can be represented 
as: 

(3.4) R(A,n) = . 

A- 

It is well know since the results established in [14, 24, 26] that the study of many properties (construction of the 
resolvent, calculation of the Fredholm determinant, construction of a system of partial differential equations for 
the determinant) of the so-called integrable integral operators I + V can be deduced from the solution of a certain 
RHP. In the case of the kernel of interest, this RHP reads 

• is analytic on C \ [ ; ^ ] and has continuous boundary values on ] ; ^ [ ; 
. xU) = o| j j jlog|^2 _ ^2| for A ^ ±q ; 

• Xi'^) - h + A~^0 (I I I uniformly in A ^ oo ; 



^1 1 ^ 

• X+U)GxiA)^X-U) for.le]-^;^[ . 
We remind that;k^+ stands for the ± boundary value ofx from the ±-side of its jump curve. 
The jump matrix (A) appearing in the formulation of this RHP reads 



(3.5) G^ (A) = /2 + 4 sin^ [nv (A)] \ ^ '^^^^ | (-e (A) , E (A)) . 



The above RHP admits a solution as long as det [I + V] 0. Indeed, it has been shown in [26] that the matrix 

■ dju . 



(3.6) x(A)=h- ■ -dfi, X~\A)^h+ ■ ' 

J H-A J fi-A 

-q -q 

solves the above RHP. This solution is in fact unique, as can be seen by standard arguments [2]. It follows [26] 
readily from (3.6) that the solution (/I) allows one to construct | F^ {A) ^ and ^ 0") |: 

(3.7) \f''{A))=x (i) \e''{A)) , ( F^ (i) I = ( F^ {A) (i) • 
3.2 Relation between and det [/ + y] 



One can express partial derivatives of det [/ + V] in respect to the various parameters entering in the definition of 
y (/(,//) in terms of the solution x {X) to the above RHP. We will derive a set of such identities below. These will 
play an important role in our analysis. 
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Proposition 3.1. Let ij > and FC^e) be a loop in U enlacing counterclockwisely '^e and such that it goes to 
infinity in the regions where e"'"^^^ ij > is decaying exponentially fast. Then, 



(3.8) 5;,logdet[/ + y] = -i—{ 

OJ] 



^ ^e'''"^'hr [(d^) (z) (cr3 + 2C [e-^] (z) o-^)x-' (z)] 



17=0+ 



(3.9) 5aologdet[/ + V] = X 



^ £ [d^uiz)] t'^^^^'hr [{d^) (z) (crs + 2C [e^^] (z) cr+)^-i (z)] 



lr(%) 



r;=0+ 



There, C [/] stands for the Cauchy transform on '^e cmd C± [G]/or its ± boundary values on '^e- One has more 
explicitly 



(3.10) C[G](^) 



J jx- X 2in ' 



and C+ [G] (A) - C_ [G] (A) ^G(A) , for A e "^e- 



Proof — The proof goes along similar lines to [32]. It is straightforward that 



(3.11) 



d, log det [/ + V] - J [d,V ■ (/ - R)] (A, A) dA . 

-1 



In order to transform (3.11) into (3.8), one should start by writing a convenient representation for dxV{A,iu). One 
has that 

(3.12) d,e (A)^-^u (A) e (A) , d,E (^) = {A) E(A)-e (A) J ^ uis) -u(A) ^_2 



2n s - A 



The last integral can be recast in a more convenient form 



(3.13) 



u{z) 



r dix u(ji)-u(A) -2. X cf) — 

J 27r n-A ^ ^' j 2n J 2in {z- A){z- n) 

.3 \ cdn r dz e'''"(^> _2, I .8 \ cd^ r 

^-%\]Yn f2rniz-AKz-,f ^ '%] i 2^: f 



e (ju) 



_ _.d_ 
drj 



77=0+ 

dz e'''"^^) rdne-^{p) 



2in {z- A) J 2n z- n 



drj 

»7=0+ 



2in {z - A){z- IJ.) 



d^e'''"(^)c[e-2](z) 



lr(%) 



27r iz-A) 



77=0+ 



We first have replaced the ratio of differences by a contour integral on r({/l,/i}). Here T ({A, fi]) consists of two 
small loops around the points A and yu. In order to manipulate convergent integrals, we then wrote the integral as 
an 77-derivative. The derivative symbol could then be taken out of the integral. Next we deformed the contour of 
integration from a compact one F ({A,/j.}) into F C^f:). Such a replacement is allowed as is chosen precisely in 
such a way so as to make e"'"^'^\ tj > 0, decay exponentially on a small neighborhood of "^e where one can draw 
F C^e). Such a choice of contours allows us to satisfy to the hypothesis of Fubini's theorem and hence permute 
the orders of integration. Also, we stress that one should compute the ?7-derivative only once that all integrals have 
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been computed. Indeed, for generic choices of functions u, permuting the 77-derivation with the /l-integration in 
the last hne of (3.13), leads to an apriori divergent integral. 

Once that this differential identity is established, one readily convinces oneself that 



(3.14) d,V(A,ix) = i— 

or] 



An 



T]=0+ 



Denoting 5 (z) - os +2C ^] (z) cr"^, using the representation (3.4) of the resolvent R and the fact that 
{F\A) I F«Ou)) = tr [I F«Ou) ) ( F^A) |], we get 



(3.15) 5;,logdet[/ + V] = /— 

or] 



d.1 



{EHA)\Siz)\EH^)) 



lr(%) 



77=0+ 



-/— tr 
dr] 



X 



1 1 



X — z X- jji 



S{z) 



\E^(M))(F^ifi)\ 



<4i-zY 



77=0+ 



Using the integral expressions (3.6) for;^^ and;^' \ we obtain 



5;,logdet[/+ V] = /|- 
077 



r {E^iA)\s {z)\e^{A)) 
' dA^ — ^ - 



r(%) 



{z-AY 



-•— tr 
dr] 
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r(%) 



Ou-z)' 



77=0+ 



(3.16) 



_ _.d_ 
dr] 



77=0+ 



where we used (3.7). The proof of identity (3.9) goes along very similar lines. 



4 The first set of transformations on tlie RHP 



We now perform several transformations on the original RHP. We first simplify the form of the oscillating functions 
E appearing in the formulation of the RHP. This step in carried in the spirit of [28]. Then, we map this new RHP 
into one whose jump matrix can be written as the identity plus some purely off-diagonal matrix. Finally, we apply 
the non-hnear steepest descent method by deforming the contour so as to obtain jump matrices that are a O (x~°°) 
uniformly away from ±q and Aq. These last steps are a standard implementation of the Deift-Zhou steepest descent 
method [16, 17]. 
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4.1 Simplification of the function E 

In order to replace the compUcated function £ by we perform the substitution 
(4. 1) X U) = X U) (h + cr^C iX)) 

where C is the rational Cauchy transform with support on '^e defined in (3.10). 

It is readily checked that this new matrix x is the unique solution to the RHP 

• ^is analytic on C \ and has continuous boundary values on ^£ \ {±^} ; 

• X^X) - I2 + A~^0 1 I I j, uniformly in /l — > cx) ; 

.^(^) = 0|J jJlog|i2-^2| for^^i^; 

• X^{A)G^{A)=x-{X) forAe^E • 
The jump matrix for Stakes two different forms 



(4.2) 

and 
(4.3) 



Gj{A) = 



for A e]-q;q[ 



Gy(A)^ 



1 e'HA) 







1 



for Ae'^eX [-q;q] 



The existence and uniqueness of solutions for the RHP for;^^ ensures that there is a one-to-one correspondence 
between and^. 

4.2 Uniformization of the jump matrices 

We now carry out the second substitution that will yield a RHP with upper or lower diagonal jump matrices whose 
diagonal is the identity. For this purpose, we define 

1 



(4.4) 



a (A) = KiA) 



A + q 
A - q 



where log 



vW)-vM 

A- fi 



d/i 



The function a (A) is holomorphic on C \ [ -q ; q], a (A) — > 1 and satisfies to the jump condition 

(4.5) a+ (A) e^''^^^'> ^a-{A) for Ae]-q;q[ . 
Then we set 

(4.6) E{A)=xiA)a''HA) . 
The matrix H (A) is the unique solution to the RHP: 



S is analytic on C \ and has continuous boundary values on "^e \ ; 
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• S(/l) = O I J \y^- qr'^'"^''^ U + qr'"^-''^ log _ for A ^ ±q ; 

• E{A) = 72 + ^"^o| j j J uniformly in /i ^ CX3 ; 

• H+(.l) Gs(A) = E_(A) for /I e ^£ . 

The new jump matrix Gn (A) reads 

/ 1 a,-2g-2 \ (10] 
Ge = ^q 1 J ^^"^EM-qiq] and = |^ ^2,vrv^^^_^2 ^g-2;;rv _ ^e]-^;^[. 

Again, there is a one-to-one correspondence between S and;^^. 
4.3 Deformation of the contour 

We now perform the third substitution that will result in a change of the shape of jump contour. Due to the fact 

that e-^{A) are exponentially small in x in appropriate regions of the complex plane, we will end up with a RHP 
for an unknown matrix T whose jump matrices are /2 + O (x~°°) and this for A uniformly away from the points ±q 
and Aq. 



4.3.1 The tune-like regime 




Figure 2. Contour St = TlT^ U r^' U rl"' U T appearing in the RHP for Y (time-like regime) . 
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We first introduce three auxiliary matrices M(A) and N^^^^\A) 
(4.7) MiA) = (J --'(Y^^^) ^I,^PiA)cr^ 



N^^\A) 



= /2 + e(^\^)(r- 



(4.8) N^^KA) 



a2(^)e2(^)(e-2'Vrva)-l)" 1 

1 ^ 

aliA) e'^''^^^'>eHA) (e-^"''^'^) - l)^ 1 



= /2 + 0^^\i)cr- . 



Note that although the matrices N^^^^\A) have different expressions, they coincide on U due to the jump conditions 

for a (A). It is clear from its very definition that N'-^\A), resp. N'-^XA), has an analytic continuation to some 
neighborhood of [ ; ^ ] in the lower, resp. upper, half-plane. Also, the matrix M (A) has an analytic continuation 
to ?7 \ [ ; ^ ] starting from ^£ \ [ ; ] . 



The functions P and Q^^^ and Q^^^ have the local parameterizations around ±q 

2mv{A) _ J 2v(A 

—7- for A e dD-q s, P(^)^ e''^' ka] 



(4.9) P (A) = e'f- J ) for A e dD-q,s , P(A)- ^"^^ "7^1^ f""" ^ ^ ^^^'-^ ' 



(4.10) 

= C(^)(i) e-'^- J'^'^ (e^--^-^) - l) and Q^^XA) = C^^\A) e'^' [^,]'''^'^ (e^''^^^'*^ - l) . 
There - x{u{A) - u {-q)) and = (m (^) - u {A}) and we have set 



[x(q-A)r^'^ [u{A)-ui-q)l '^ 



We now define a piecewise analytic matrix T according to Fig. 2. We will be more specific about the choice 
of the contours T^^^^^ around the points ±q and Aq when we will be constructing the local parametrices. Here, 
we only precise that the jump contour for Y remains in U and that all jump curves are chosen so that, for a fixed 
z e rf^ U rf ^ \ {±q,Ao] (resp. z e rf^ U r|^^), e'-'"^^^ (resp. e-'^"^^)) is exponentially small in x. The matrix T is 

discontinuous across the curve = F^^^ U F^^^ U F^^^ U F^^\ One readily checks that the matrix T is the unique 
solution of the below RHP (and hence there is a one-to-one correspondence between and Y): 

• Y is analytic on C \ and has continuous boundary values on 'Lr \ {±^} ; 

• Y(i) = o|j j J (^ - ^)-'^3y(?) + ^)(r3y(-?) log _ A^±q; 

• r(A) + A-^O I j j J uniformly in /I — > oo ; 

• Y+(^) GriA) = Y_(^) for ^ e St \ {±q, Aq] . 

With 
(4.12) 

Gy{A)^M{A) on fJ,^^ U fJ^^ , Gy{A)^N^^\A) on vf and Gy {A) ^ N^'^Xa) ontfK 
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Figure 3. Contour 2t = ^ U T^' U r!" U r^' appearing in the RHP for Y (space-Uke regime). 



4.3.2 The space-like regime 

We introduce two matrices M and 



1 a e 



2^-2 \ ( I ^ 



(4.13) M = ( Q " ; j^l2 + P(A)<r- and N = \^ _ ^ j 



= /2 + e(/l)cr- . 



The matrix M(A) has an analytic continuation to ?7 \ [ ; ^ ] starting from \ [ ; ^ ]. The matrix N(A) has an 
analytic continuation to U nB.-. 

This allows to write convenient local parameterizations around ±q for P and Q: 
(4.14) 

P(.) ^ e'-^-. k-J-"^^^ -"t' " ' - -e'-^- k/''^ '"T " ' with j^- - "^?]r"^7f? . 

c('f-)(^) L^^J CW(i) I = x(m(^)-m(^)) 

Similarly, 

(4.15) e(^) = C(^)(A) e-'f- [^_J''^'^(e2-^W) - l) = -C(«>(A) e-'^^ [^J"'"^"^ (e-'''-^^^ - l) . 



Here, we bring to the reader's attention the difference of signs in the definition of in the time-like and space- 
like regimes. Also, the functions C'-^^'^XX) have been defined in (4.1 1). The sole difference is that, in the space-like 
regime, the +iO^ regularization plays no role. 

The matrix T (A) defined in Fig. 3 in the unique solution to exactly the same RHP as formulated for the time-Uke 
case but with the contours being defined in Fig. 3 and the jump matrix being now given by 

(4.16) Gr(A)^M(A) on F^^^ U fJ,'^^ U fj^^ and GriA)^N(A) on F^^ 
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5 The local parametrices. 

We now build the parametrices around ±q and Aq. These will allows us to put the RHP for T (and hence the 
one forx) in correspondence with a RHP that has its jump matrices close to the identity, uniformly on its whole 
jump contour (in the case of T (A) the jump matrices are close to the identity only uniformly away from the points 
Aq and ±q). The role of the parametrices is to mimic the complicated local behavior of the solution x near the 
stationary point Aq and the endpoints +q. Once again, due to slight differences between the two regimes, we treat 
the space-like and the time-like regimes separately. 

5.1 The time-like regime 

We recall that for the time-like regime, the functions P (A) and Q^^^^\A) appearing in the jump matrices are given 
respectively by (4.9) and (4.10) with C'^^^'^XA) given by (4.11). 

5.1.1 The parametrix around Aq 

It follows from the assumptions gathered in subsection 2.1, that the function u admits a local parameterization 
around Aq, ie there exists 6 > Q such that i^Ao,s c U and a holomorphic function h on some open neighborhood of 
DAf,,6 such that u{A)-u {Aq) = -o? {A) with oj {A) -iA-Ao)h (A), and h (Dao,6 n H±) c H±. 

The curves ^f'/f^ in i)^,i5 are defined according to Fig. 4. 




Figure 4. Contours appearing in the local RHP around Aq in the time-like case. 

The parametrix Pq around Ao reads 

(5.1) ^o(^) ./2-^2i(A)e-'t^^^^|^^(l,^;ixa>2(A)) a- . 

There, *P (a, b; z) is Tricomi's confluent hypergeometric function whose definition is recalled in appendix A. 
The function b2i is defined piecewise: 

(5.2) b2i(A) = a2(i)e-'^"^'^«>-«^'^^(e-2''^^^'^>-lf for A e BT n , 

(5.3) b2l(A) = ^2(^)g4/^(a)g-,>«(^)-gU)Jg-2;;rva)_ fo^ /I £ H ©^,5 . 

It is holomorphic on £)io,5 due to the jump condition satisfied by a (4.5). 
The paramertix Pq solves the RHP: 
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• Po is analytic in Dao,s \ {^f ^ U r[^^} n Dao,6 with continuous boundary values on |r^^^ U T^^^} n D^^^s ; 

• ^0-12 + —pO (o-~) uniformly in A e dD^^s ; 

The first two points in the formulation of the RHP for Pq are obvious. The vaUdity of the jump conditions can 
be checked with the help of identity (A.l). 

5.1.2 The parametrix at -q 

The parametrices for the local RHPs at ±q are well known. They have already appeared in a series of works 
[4, 11, 32] and can be constructed from the differential equation method [22]. Here, we recall their form. 

The parametrix P-q around -q reads 

* (y (i) , 1 ; -ix [u (A) - u {-q)]) ibniA) {\ - v{A) ,l;ix[u{A) - u {-q)])\ 



(5.4) P_q (A) = T (/I) L (A) [x (u (A) - u {-q))Y^^^'^' e" 



(5.5) ^{A) = 



-ib2\{A) '^{\+v{A),\;-ix [u (A) - u {-q)]) Y (-v(^) ,l;ix [u (A) - u {-q)]) 



.sin[;Tv(i)] .... 
(5.6) 1rV)(i) ' ^^t^^t bn{A)b2i{A) = -v\A) . 

^^^^^^ ~ ~^in[;ry(^)]r2(-y(^)) 
C'^^KA) is given by (4.1 1) and 



(5.7) L{A) = 



1 

e- 



-njl < arg [u{X) - u{-q)\ < njl, 
i^(X) j < arg [u {A) - u {-q)] < n, 

-7T < arg [u (A) - u (-q)] < -njl. 



g2;^(,l) Q 

1 



P-q is an exact solution of the RHP: 



• P-q is analytic on D-q^s' \ {r^^^ U F^^^j with continuous boundary values on |r^^^ U F^^^j \ {-q} ; 
. P-q{A) = O I J j j (^ + qY'^^-i^ log \A + q\, A^-q; 

• P-q {A) -I2 + 1 I I j ' uniformly in A e dD-q^s' \ 

[p-q]^{A)MiA) = [p-q]_{A) iorAeYf^nD-q,s'; 
[p-q]j,A)N^^\A) = [P-q]_W ioxAeTfr\D-q,5' . 
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Here, we have set 

(5.8) = 2sup {\% (v(A)) : A e dD^g,3' U dD-g,s'\} < 1 . 

The fact that ps' < 1 for sufficiently small 6' is a consequence of the assumptions that \% (v (±^))| < 1/2. The 
canonically oriented contour dD-q^s' together with the definition of the contours r^^| is depicted in Fig. 5. 5' >Q 

o 

is chosen in such a way that D±q^5' c U, D±q^s' ^ Oao,s = and Dq^g' n D-q^s' - 0. Playing with the 6 entering 
in the definition of the parametrix Pq, one can tune it in such a way that 6' - 5. We shall assume such a choice in 
the following. 

u - u (-q) 




dD-q,6' 



Figure 5. Contours for the local RHP around -q in the time-like case. 



5.1.3 The parametrix at q 

The parametrix Pq around q reads 



(5.9) 
Here, 

(5.10) T(A) = 



Pq {A) = Y (i) L (i) [x {u {q) - u {A))] e" ^ 



'¥{-v{A),\; -ix [u {A) - u {q)]) 
-ib2\ {A)'^{\-v{A),\; -ix [u (A) - u (q)]) 



i bn(A) '¥{\ + v(A),\;ix[u(A)-u (q)]) 
'^(v(A),\;ix[u(A)-u(q)]) 



(5.11) 



bniA) = 



^[C(«)(i)]"' 



r2(-y(^)) sin [Try (i)] 
hi (A) - OT-ir2 (1 - V (A)) C(^)(^) sin [nv (A)] 
C^'^\A) is given by (4.11) and 



, bn{A)b2i(A)^-AA) . 



(5.12) 



L(A)^\ 



h -nil < arg [u (q) - u (A)] < njl , 

1 



1 

e- 



n/2 < arg [u (q) - u {A)\ < n , 
2^(X) j -n<dxg [u (q) - u (A)] < -n/l 



Pq is an exact solution of the RHP: 



Pq is analytic on£)^_5\|r^^^ur^^^|n£)^_5 and has continuous boundary values on |r^^^ur^^^|n£)^ ,5\{^} ; 
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• "PqU) = h + ^1^0 ^ J J j , uniformly in e dDq^s ; 
' [rg]^(A) N^''\A) = [rg]_(A) for A e rj^^ n \ {q} , 
M{A) = [Pg]_{A) for /I e rf ^ n D,,^ \ iq) ■ 



The canonically oriented contour dDq^s together with the definition of the curves r^^| in the vicinity of q is 
depicted in Fig. 6. Note the change of orientation of the jump curve due to u' (q) < 0. Also p<j is as given in (5.8). 



u — u{q) 




Figure 6. Contours for the local RHP around q in the time-hke case. 



5.1.4 Asymptotically analysable RHP for n 

We now define a piecewise analytic matrix IT in terms of Y and the parametrices according to Fig. 7. In particular 
one has IT = T everywhere outside of the disks. The matrix IT has its jump matrices uniformly close to the identity 
matrix in respect to the x +oa limit. Hence, it can be computed perturbatively in x by the use [17] of Neumann 
series expansion for the solution of the singular integral equation equivalent to the RHP for H. This matrix H is 
the unique solution to the RHP: 

• n is analytic on C \ Xn and has continuous boundary values on En ; 

• n(^) = /2 + A-^0 1 j j j , uniformly in /I — > oo ; 

• n+(^) GnU) = n_(i) for ^ e En . 



The jump matrix Gn {A) for H {A) reads 



(5.13) Gn {A) = Gt {A) on F = rf ^ U F^^^ U vf^ U F^'^^ and Gn {A) - 



rilU) on 



P-'{A) on -dDA,,5 



5.2 Asymptotic expansion for the algebraically small jump matrices 

Note that the jump matrices along F are exponentially close to I2 in x and this in the L} (En) n (En) n L°° (En) 
sense. Only the jump matrices on the disks are algebraically in x close to the identity matrix. The latter jump 
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Figure 7. Contour Sn appearing in the RHP for n, time-like regime. 



matrices have the below asymptotic expansion into inverse powers of x, valid uniformly on the boundary of their 
respective domains of definition (dD±q^s or OD^^s)'- 

V(-'"\s) 



^(n + l)![x(5 + ^)ri 



(5.14) ^o'(^) - /2 + (^-y — 



Where, 



\u{s)-u{-q)j \ 1)^^21 (1 +v)„ (y)„_^i / 

We remind that o) {A) - {A - Aq) h (i) and we have used the conditions det \P±q\ = 1 = det YPq\ so as to invert 
the parametrices and then infer their asymptotic expansion from the one of CHF (A.2). Also, we have not made 
explicit that bij, bjj and v are functions of s. 



5.3 The space-like regime 



The discussion of the space-like regime resembles, up to minor subtelties, to the previous one. Therefore, we 
make it as short as possible. 
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5.3.1 The parametrix around Aq 

The parametrix on Da(,,6 for the local RHP around Aq reads 

(5.18) PQ(A)^l2-bnU)e''^^!^^^^'¥(l,l--ixoj\A)]o-^ with bn(A) ^ a-\A) e''"^^^*'^^\ 

2in \ 2 / 

^712 is holomorphic on 0^0,6- The parametrix is a solution to the RHP 

• Pq is analytic in Da^^s \ '^f^ U T^^^j n D^o,^ and has continuous boundary values on jp^^^ U P^^^j n Daq,6', 

• Pq- h + —p^ (c*) uniformly in /I e dD^^s ; 

The jump curve for the parametrix Pq is depicted on Fig. 8. 




Figure 8. Contours for the local RHP around Aq in the space-like case. 



5.3.2 The parametrix around -q 

This parametrix P-q is exactly the same as in the time-like regime. Hence, we do not present it here. 

5.3.3 The parametrix around q 

The parametrix Pq around q reads 

(5.19) Pq {A) = ^ (^) L {A) [x {u {A) - u iqW^^^'^' e^ . 

Here, 

'^{-v{A),l\-ix[u{A)-u{q)]) /^i2 (^) 'P(l v(/l) , 1; « [m (i) - m (^)]) ^ 

-ibiiiA) {\ - v{A) ,\;-ix[u{A) - uiq)]) {v{A) ,\;ix[u{A) - u{q)\) 



(5.20) W{A) 



W n(K)(i\ sm[;ry(/l)] _ _ 

(5.21) _ inC^^\A) ' bn{A)b2M) = -v\A) 

~ P2(y(^))sin[7ry(^)] 
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C^R) is given by (4.11) and 



(5.22) 



L{A) = 



h -nil < arg [u {A) - u {q)] < njl , 

I Q-ihzvU) J < arg [u (A) - u (q)] < n , 



g-2OTV(/l) Q 

1 



-n < arg [u {A) - u (q)] < -7t/2 



fa is an exact solution of the RHP: 



• Pq is analytic on Dq^s \ U F^^^} n Dq^s with continuous boundary values on |r^^^ U rj^^} n Dq^s \ {^} ; 
. rq{A)^oi^ j j j(^-^)-'^3v(?)iog|^_^| ^ x^q; 

• "Pqi^) -h + ^I^O I J J j , uniformly in .1 e dDq^s ; 

' [pq]^{X) M {A) = for A e if^ n !Dq,s , 

• I 

[pj^(^) {A) = [nlj-l) for A e T^^^ n Dq,6 . 

The canonically oriented contour dDq^s as well as the definition of the jump curves T^^^^^ is depicted in Fig. 9. 
Finally, ps has been defined in (5.8). 



u- u{q) 








Figure 9. Contours for the parametrix around q in the space-hke regime. 



5.3.4 The RHP for H 

The matrix 11 is defined according to Fig. 10 and is the unique solution to the RHP formulated in exactly the same 
way as for the time-Uke regime. The difi'erence consists in the precise form of the contours due to the fact that in 
the space-hke regime Aq > q. 



5.4 The asymptotic expansion for the parametrices 

The jump matrices !P+^ have the same asymptotic expansion as in the time-like regime (5.14) with the sole excep- 
tion that the coefliicients bi2, ^21 entering in the definition of V^*''"^ (5.16) are now given by (5.21). The matrix 



24 



RIEMANN-HILBERT APPROACH TO THE TIME-DEPENDENT GENERALIZED SINE KERNEL. 




Figure 10. Contour Sn appearing in the RHP for IT, space-like regime, 
f Q ^ has the below asymptotic expansion 

(5.23) (^) - - Z : ''"'^'^ . with .(«)(.) = -tnis) (-if '^^'"'^J:,, 

^ y'+J (5 - ^)2"+i 2;r /i2n+i(j) 



6 Asymptotic expansion of the Fredholm determinant 



Starting from now on, we will treat both regimes (space and time-like) simultaneously. 



6.1 The asymptotic expansion for IT 



In this subsection we present two ways of writing down the asymptotic expansion for the matrix IT. The first, 
given in proposition 6.1, traces back all the different fractional powers of x and oscillating terms that appear in 
the asymptotic expansion of IT. It also provides one with a sharp and quite optimal control of the remainders. The 
second one, given in proposition 6.2, is considerably less exphcit and, by far, does not provide optimal estimates 
for the remainders. However, it is easier to implement form the computational point of view, especially when one 
is interested in calculating only a couple of terms in the asymptotics. One can then build on the first asymptotic 
expansion so as to on the one hand argue for a sharper form of the estimates for the remainders and on the other 
hand identify which among the computed terms are relevent and which are not. We start this section by presenting 
the Neumann series expansion for H. 

Definition 6.1. Let Sn be the jump contour for the matrix H. We define the contour as being the inslotted 
version of the N-fold Carthesian product Sn X • • • X ^n- Namely it is obtained from Sn X • • • X En putting the 
contour for zk+i slightly shifted to the right from the contour for Zk- We have depicted the inslotted contour for 
N = 2 on Fig. 11. 

Let pr^. stands for the projection on the k^^ factor of an N-fold Carthesian product, ie given z = (zi, ■ . ■ ,Zn) 
one has pr^ (z) = Zk- The contour thus defines N curves Xn {Zk\ = Wk (^n'^)' ^ ~ 1» • • • > Each of these 
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Figure 11. The inslotted contour for = 2. The integration over z.\ runs through the dotted 
contour whereas the one over zi runs through the full one. dO {zi\ refers to the three disks over 
which the variable Zi is integrated. 



can be interpreted as giving rise to the jump contour for the RHP problem associaed with the matrix IT. In the 

following whenever A, resp. V, is inegrated along En {_Zk\, it should be understood as originating from the jump 
matrix h + A, resp. I2 + V, appearing in the RHP for IT, resp. n~\ when the latter is formulated on the jump 
contour En [zk]- 

Lemma 6.1. Let I + Abe the jump matrix for H and V = Comat (A)' be the transpose of the adjugate matrix to A. 
Then, for x-large enough, the matrices H and admit the below uniformly convergent Neumann series 



(6.1) n(^) = /2 + 



(6.2) U-\A) = I2 + 



d^z A(z^)...A(zi) 

A^>1^, ^^^^^ ^X-z,) n iZs-Zs^l) 



2J 

A'>1 «n 



d^z V(zi)...V(Zjv) 

. {2in)\^ /-I , 
m (/l-zi) n iZs-Zs+\) 

s=l 



The convergence holds in U° (O) sense for A e O, with O any subset ofC such that d{0, En) > 0. Also, it holds 
for A± € En in the L^(En) sense. Finally, the matrices A and V that are integrated along the inslotted contour e|^^ 
should be understood according to definition 6.1. 



Proof — 

We define two linear operators on 2 x 2 L^(En)-valued matrices 

/ds r ds 
— — -M is) A is) and *Cl^ [M] {X) - — — V {s) M {s) . 
lin{A+ - s) "J lin{A+ - s) 

Using that for sufficiently regular, not necessarily bounded, contours En, the ± limits of the Cauchy transform 
with support on En are continuous operators on L^(En) with norm c (En) [31], it is easy to see that that the two 
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above operators are also continuous-^ on the space ^2 (^■^(^n)) of 2 x 2 matrices with L-^(£n) entries: 
(6-4) ||ct^ [M]|| 

(6.5) ll'Czn M|L2(Sn) ^ 2c(2:n)||A||z.~(Sn)l|M|L2(s„) 

There we made use of the fact that V is the transpose of the adjugate matrix to A so that HAH^oo^j^^^) = ||V||^co(j;jj) 
and ||A||i2(2n) = IIV|lL2(Sn). 

It is a standard fact [5] that there is a one-to-one correspondence between the solution to the RHP for n (or 
n~^) and the unique solution to the singular integral equations 

(6.6) n+-C^jn+]=/2 and Cl^[u-'] ^ h ■ 

Indeed, provided that 11+ is known, the matrix 11 (or n~^) admits the below integral representation for A away 
from Sn 

(6.7) U{A) = l2+ f ^. n+ is) A is) and U-\A) = h + f ^. V (s) U^' (s) . 

J 2m {A - s) J 2m {A - s) 

The estimates for the jump matrices on the boundary of the discs dO±q^s and dOxQ,5 and the specific choice for 

the shape of the contour En at infinity lead to ||A||^2(j;jj) + l|A||/,oo(j;j,) - O (x"*) with w = min (1 /2, 1 - ps) > 1 and 
defined in (5.8). This implies that for x-large enough the operators / - and / C^^ are invertible and that 
their inverse can be computed by a Neumann series expansion converging in (Ln)'- 

(6.8) n.(.) . /2 . 2 Kn^^] u) = ^ z / £^ "'t[-"''' • 

(6.9) n-(.) - / £^ '^tv'^"^ • 

s=l 

Where {C^jj}^= C^^ o ■•• o C^^ stands for the composition of operators C^^. In (6.8)-(6.9) the integration runs 
across the Carthesian product of copies of En: Sn x • • • x Sn. 

The fact that 11*^ (A) admits a uniformly convergent Neumann series for A belonging to any open set O 
at finite distance from En follows from the (Sn) convergence of the series (6.8)-(6.9), the fact that A e 
^2 (l^ C^n) n (En)), and that d (O, En) > 0. 

Finally, it is easy to check that one gets the expression for 11(A) (resp. 11"^-^)) on C \ En by replacing the + 
type regularization A+ of A in (6.8) (resp. (6.9)) by A e C \ En. 

The A^"^ summand of the Neumann series for 11*' can be expressed in terms of regularized by deforming the 
original contour En X • • • X En to the inslotted one e|^\ The latter manipulation is possible due to the properties 
of the locally analytic matrices A (z) and V (z). It allows one to get rid of the + regularization in the integrals. 

The construction of the inslotted contour e|^^ is depicted in Fig. 12 and 13. Initially, the integral is performed 
with the use of the + boundary value of zi on the integration contour for zi- Hence, away from the points of triple 
intersections c,, we can deform the integration contour for zi to the + side of the integration contour for Z2- One 



^By intechanging the roles of A and M, it is easy to see that C^^ is continuous on ^2 {L°°(^n)) since A 6 L^(En). 



K.K. KOZLOWSKI 



27 



ends up with a contour as depicted in Fig. 12. There, the dotted Unes correspond to the integration contour for zi 
whereas the full lines give the integration contour for zi- One then proceeds inductively in this way up to zat. As 
A is assumed to lie uniformly away from the original contour En, there is no problem to deform the integration 
contour for zi in the vicinity of En as the pole at zi = A is lying "far" away. 

It remains to threat the integration on the intersection points of the disks 5D[zy] with the curves rj,^^^^[zj]. 
We first reduce the most interior disc (the one over which zn is integrated and then the procedure is repeated by 
induction) to smaller a one. The jump matrices A have different analytic continuations from the right and left of the 
points Ci (this corresponds to the discontinuity lines of Pq and Taking this difference into account produces 
the small extensions of the contours T^^^^^ [zn] as depicted on the right part of Fig. 13 together with smaller discs 
dD [zn]- It is in this way that the matrix A integrated over En [zn] is identified with the one stemming from the 
jump matrix for H when the latter is defined as in (5.13) but with jumps on En [zn] (what corresponds to slight 
deformations of the curves rl^(^^). ■ 




Figure 12. Construction of the inslotted contour. 





Figure 13. Deformation of the circles. 
Proposition 6.1. The matrix li admits the series expansion 

(6.10) - - ^ -^"-^^^ 



n(i) = /2 + 2 



N>1 
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that is valid uniformly away from 2n and also on the boundary Sn in the sense of I?' (Sn) boundary values. The 
coefficients Xl^ of this expansion take the form 

[W2] ixm[u(q)-u(-q)\ 

(6.11) W)^A,iA)^ ,24v(.).v(-.)] 

[Nil] h i'^yb i^,n[u{q)-u(-q)\ ixt,[bu{Ao)^pu{q)+{p-h)u(-q)\ ^ ^ ^ 

+ V V V - xi- n^'"'^'^V/i) 

Z-l Z-l Z-l Jlm\v{q)+v{-q)\ Jlt,{b-p)v(-q)-lpijv(q) ^^A' ' 

b=l p=0 m=b-[§] ^ 

and one should set i] - I in the space-like regime and t] — -\ in the time-like. 



The matrix An{A) contains only exponentially small corrections, ie [Ajq\j{X) — 0{x °°) with a O that is 
uniform for A-uniformly away from Sn- 



The matrices n^\/l) and n^'*'^\/l) admit the asymptotic expansion 



(6.12) 



nr(i) = X^^Z^X) with nj!;(^) = o 



r>0 



i^ogx) 



N+r-6m!3-'2-m 



-M 



(6.13) nf'''P\A) = J]nj!;^'^\i) with n^Zr^p\A)^o 

r>0 

The estimates hold for A uniformly away from £n- 

The matrix M appearing in the various O estimates takes the form: 

^ixuiq) 



(log 4 



N+r-2(m+b) \ 

M 



(6.14) 



M: 



1 



m 



m. 



2v{q) 



_ e-'^"(-«' _ ^/x 



+ x2y(q) 



ixu{Ao) 



+ m_ 



r2v(-4') 



+ mo 



g;;cM(/lo) 



1 



There m±, mo, m± and mo are x-independent coefficents. Moreover, necessarily, mo - in the time-like regime 
and ifio -0 in the space-like one. 



We postpone the proof of this asymptotic expansion to appendix B as it is rather cumbersome and long. How- 
ever, at this point, we would like to make several comments in respect to the form of the asymptotic expansion. 

The above asymptotic expansion is in a form very similar to the one of the functionals Hn [v, e^, u] given in 
Theorem 2.2. In fact, the large-x behavior of the matrix Il^f contains various fractional powers of x, each appearing 
with its own oscillating pre-factor. Once that one has fixed a given phase factor and fractional power of x, then 
the corresponding matrix coefficients nj^"^ or nj^'' ^' admit an asymptotic expansion in the more-or-less standard 
sense. That is to say, each of their entries admits an asymptotic expansion into a series whose n* term can be 
written as Pn+u (log.^:) /x" with Pn+h being a polynomial of degree at most N + n. One of the consequences 
of such a structure is that an oscillating term present in x~"n„ (A) may be in fact dominant in respect to, say, a 
non-oscillating term present in x~" n„' (A) where n' < n. 

We would also like to point out that the asymptotic expansions of n^v and n^v+i share many oscillating terms at 
equal frequencies (eg Q'"["('i'>-"(-y)] is present in Hn and Hn+i for any > 2). However, those issued from IIn+\ 
have an additional dumping pre-factor log;ic ■ x~^ in respect to the same ones issued from H^r. Finally, there may 
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also appear additional oscillatory terms e*'-^"^^\ z = ±^ or Aq (and their associated fractional powers of x) in the 
off-diagonal parts of nj^ and nj^'*'''^ cf (6.14). 

There is also another way of writing down the asymptotic expansion of n (A). Although it is more compact, it 
is also less explicit and provides one with weaker estimates for the remainders. 

Proposition 6.2. The matrix Yi admits the asymptotic expansion 

(6.15) n(^) = /2 + |;n!^+o( j j with w = min(i I-P.^) , 

„>() X ^ \ I \ I 

that is valid uniformly away from En- 

For A belonging to any connected component of oo in C \ Sn, the first few terms appearing in this expansion 
read 

(6.16, = n<..a, = -2^, 

A-Aq A-eq 



(617) U^^HA) = y d^'^KAo) i V^'''>\eq)a (rV(^'"W ) cr 

Aq- eq \ A - Aq A — eq J 2 ds^ \ A — s ] 

The expression for 11^^^ is a bit more involved. 



1 V — ^ ^^''^^ ~ 2V(^'*'\e^) (j) 



2^ds[ A-s , ^^^^ 

There a = a'^ in the space-like regime and cr - cr~ in the time-like regime. 

This form of the asymptotic expansion is the closest, in spirit, to the one appearing in the literature, cf eg [15]. 

However, it does not represent a "well-ordered" asymptotic expansion in the sense that each matrix 11^"^ depends 
on the various fractional powers of x and oscillating corrections. Some terms present in the entries of n*^^'^ are 
dominant in respect to the ones present Yl^^\ { < p. Moreover, the expansion (6.15) does not provide one with a 
precise identification of these terms. This form is however very convenient from the computational point of view, 
and having explicit expressions for the matrices Yl^^^ easily allows one to identify the various matrices entering in 
the "well-ordered" asymptotic expansion (6.10)-(6.11). 

Proof — 

The unique solution 11+ to the singular integral equation (6.6) equivalent to the uniquely solvable RHP for H 
provides an integral representation for H on C \ Sn. Namely, 

(6.19) UiA)=l2+ f— ^— n+(^)A(^) for ^ e C \ In • 

J 2m (A - s) 

The only places where the jump matrix for H is not exponentially close to the identity are the three boundaries 
of the discs -dD±q^s and -dI)Ao,6- There one has 

(6.20) ^^"Pll = n_ on - dD±q,6 and n+pQ^ = n_ on - dD^^^s ■ 
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Note that the minus sign refers to the clockwise orientation of the boundary of the discs in Fig. 10 and 7. 

By using the estimate 

(6.21) W(A) = ||A||ii(Sn) + ||A||i»(Sn) = 0(x-"') with w ^ mm^^,l - ps"^ 

and equation (6.6), one shows by standard methods (see eg [15]) the existence of the asymptotic expansion (6.15) 
for n (A). This expansion is valid uniformly away from the jump curve Xn. 

We would like to stress that for computing the coefficients of the asymptotic expansion, we can drop the 
integration contours other then the boundaries of the disks dD±q/;ig;s. Indeed as (11+ - 12) e (2n), cf (6.1) , 
and l|A||^2p^i(^f^ - 0(x '"), it is clear that the integration along T in (6.6) can only produce exponentially small 
corrections. Thence, it cannot contribute to the asymptotic expansion (6.15). As a consequence, the matrix 
coefficients n*^"^ in (6.15) can be computed by plugging^ the asymptotic series into the integral equation (6.19), 
dropping there all the exponentially small corrections (stemming from the integration along T) and replacing ffie 
jump matrices !P~^, !P+^ by their asymptotic expansions which are valid uniformly on the boundaries of ffie three 
discs. This leads to the formal (in ffie sense that valid order by order in x) equation. 



(6.22) 



u(A)-i - — f y n+(^)^^"^(^)^ - — y r y _ ^^^'^ ^^''"^^^^ 



It now remains to equate the coefficients of equal inverse powers in x. This yields sets of recurrence relations 
between the various terms appearing in the asymptotic expansion for 11. A straightforward residue computation 
leads to ffie result for n^''\ n = 0, ... ,3, for A belonging to any connected component of cx5 in C \ En. ■ 



6.2 Proof of the leading asymptotics of the determinant 

We now prove ffieorem 2.1. We divde the proof into ffiree part. First, we obtain a modified version of ffie integral 
representation (3.8) for dx log det [/ + V] that will be more suited for our further computations. Then, we use this 
integral representation so as to compute the first few x-dependent terms in ffie asymptotics. Finally, we fix the 
constant, x-independent part of the asymtptotics. 

• Modification of the integral representation 

The first few terms of the asymptotic expansion of det [/ + V] can be obtained by using ffie identity (3.8) between 
the x-derivative of log det [/ + V] and ffie RHP datax, together wiffi ffie asymptotic expansion for 11. As a starting 
remark, we observe that one can always choose the contour TC^e) appearing in (3.8) in such a way ffiat it only 
passes in ffie region where 

(6.23) XW-^ a''"' iX) [h + C [e'^] {A) cr+) . 

Then, plugging this exact expression for^ into the trace appearing in (3.8), one gets ffiat 

(6.24) ti- [d;a iX) [cr3 + 2C [e-^] {A) (r^]x-\X)] = tr [S^H {A) a^nr^ {A)] - 2^;^ (log a) (A) . 

It remarkable, but also important from the computational point of view, that the matrix allowing one to simplify 
the comphcated functions E (A) (1.2) appearing in the formulation of the initial RHP, does not play a direct a role 

'It is possible to insert the asymptotic expansion, whicii a priori is valid only uniformly away from Sn in (6.15) in as much as one 
slightly deforms the contour En in the + direction what is allowed in virtue of the analytic properties of n+ and the local analyticity of A. 
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in the computation of the asymptotics of the determinant. In particular, one does not have to deal with integra- 
tions on T{^e) of Cauchy transforms C {X). Inserting (6.24) into (3.8), one obtains that the contribution of 
-2dx (log a) {X) can be separated from the rest, so that 



(6.25) log det [/ + V] [v, w, g] - a_i - i— 

OJ] 



J^e'''"(-^)tt[5an(^)cr3n-'(i)] 



lr(%) 



where we have set 

h;i l(y (^\\ r 

-q -q 



1 1 

(6.26) a_i = J ^ u' {A) log J - / J u'{X) v ( 



Note that one cannot exchange the //-derivation and the /l-integration symbols in (6.25) yet. To be able to do 
so, we deform the most exterior parts of the contour T{^e) in (6.25) to Y^^^ and Y^^, cf Fig. 14. We denote ■)^^^ 
the resulting interior loop. The integrand along y^*^^ remains unchanged. However, when integrating along the 
contours one should replace tr [d^II (/I) 0-31!"^ (/l)] by the difference between the two boundary values: 

(6.27) ti- [5an_ {X) a-^TT} {X)] - tr [5^n+ {X) crgH^i {X)] . 

We remind that not only the ± boundary values themselves, but also theses of II's derivatives do exist on Y^y^^- 

This is a consequence of the fact the jump matrix /2 + A for 11 on Y^y^^ admits an analytic continuation to a 

neighborhood of these curves. This fact allows one for a local deformation of the jump contour Zn, meaning that 
n+ (resp. n_) admits an analytic continuation to some neighborhood of En located on its - (resp. +) side. 

The difference in (6.27) can be estimated with the help of the jump condition for 11 along Y^y^^'- H+M = n_ 
where M is given by (4.7) : 

(6.28) ti- [5^n_ {X) a-^Xr} {X)] = tr [{[5^n+ {X)\ M {X) + n+ {X) dxM {X)} o-^M-^ (X) U-^ (X)] 

= tr [n^i (5^n+ (X)) (X) 0-3] + ti- [d;iM (X) o-^M-^ (X)] . 
Using that M - I2 + Pcr^, with P being defined in (4.9) , we obtain the jump formula 

(6.29) tr [5an_ {X) cr^ir} {X)] - tt [^;^n^ {X) cr^Yll^ {X)] = 2a-^ {X) e-^{X) tt [a^n^ {X) cr+n;i {X)] . 

Using, once again, the jump condition on Y^y^\ we see that tr [54n(/l)cr'^n"^ (/l)j has no discontinuity across 

those parts of r^^^^^ that we focus on. It can thus be extended to a holomorphic function in some neighborhood 

of this curve. Thence, we can deform the contours of integration Y^y^^ to y^^/^^ as depicted on Fig. 14. Once 
that this has been done, there is no problem anymore to exchange the 77-derivation with the /l-integration. Indeed, 
tr [5^n (/I) cr"^n (/I)] is bounded when % {X) — > ±cx5 along y^^l^\ and the function G given by 

(6.30) G{X)^u U) {l^(o)(A) + 2a-2 {X) e'^ {X) V)uyw('l)} 

is integrable. Here, \a stands for the characteristic function of the set A. Once that the 77-derivative is computed, 
we get the below integral representation 

(6.3 1) log det [/ + y] [v, u,g\^a-i+j)^G {X) tt [S^H {X) a {X) XT^ (^)] . 

y 

The final contour 7 is depicted in Fig. 14 and the matrix-valued function reads cr {X) - cr3ly(0)(/l) + cr"^l^(L)uy«)(/l). 
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Figure 14. Contour 7 = 7^^^ U 7*^'^^ U 7^*^ The contour Xn is depicted in dotted lines. 



Extracting tlie first few x-dependent terms 



tr[5anU)cr+n~H^)] is bounded on 7^^^ U 7^^^ and ||G||i,i(yj?)uyL)) = 0(x~°°). Hence, we can drop the part of 
integration over 7^^^ U 7^^^ when computing the asymptotic expansion of logdet [/ + V]. It thus remains to treat 
the integration along 7^*^\ 

As follows from proposition 6.1, IT has a uniform asymptotic expansion on 7'^*'^ given by (6.10). In order to 
obtain the leading asymptotic expansion for the .^-derivative of the determinant, it is readily seen that it is enough 
to plug in the more compact expansion (6.15) to the desired order and then drop all the terms that are irrelevant. 
This is simpler from the point of view of computations and justified a posteriori by the form of the well ordered 
asymptotic expansion (6.10). Therefore we get 

(6.32) tr [d^U (A) o-sU-' (A)] = -tr {[u^'^]\a) era) + ^tr [[n^^^]\A) - U^^\A) [u^'^]{A) - n^'\A) [n(0)]'(i)) ^3 



X2 



+ 



ltr{[rf3)]'(i)-n(>)(^)[n(i)]'(^))cr3 + o 



^±ix[u{q)-u(-q)] Qii]x[u(Aa)-u(±q}] ^^Qg x) ' 



^±2[v(?)+2v(-9)]+2 ' ^|t2v(±9) ' 



uniformly on 7^*^^ There the o refers to sub-leading terms that have been ignored. It distinguishes between the 
various oscillating and non-oscillating corrections that have been ignored. Also one should set 7] = 1 in the 
space-hke regime and tj - -1 in the time-like regime. Note that due to the compactness of 7*^°^ the order of the 
o-remainder is preserved by the integration along 7^"\ 

Note that, in (6.32), we have been able to simplify certain products by exploiting that, regardless of the time or 
space hke regimes [iT^*^^] = and that traces of matrices proportional to cr (with o" = cr* depending on the space 



ot time-hke regime) vanish ( eg H'^^^II^^^ |^n*^*^^J oc cr). 



We now insert the explicit form of the first few matrix coefficients appearing in the expansion of IT and then 
integrate the expansion (6.32) along 7*^^^ with the appropriate weight. At the end of the day, by using the precise 
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estimates provided by the expansion (6.10), we get 

osc / /, „\\ ^no / /I W / „ \ ' a'^^'^ 



(6.34) ai . i y d^o\Ao) ''^f'^ ' ""['f tr{v^-'\eq) [<r,,cr]} 



Above the last O corresponds to higher order oscillating correction with bigger phases than those involved in the 
definition of ai and a°^^. The term responsible for the logarithmic contribution to the determinant coincides with 
the one appearing in the time-independent case (the so-called generalized sine kernel) considered in [32]: 

(6.33) ao^-(v\q) + y\-q)) . 

The first Aq dependent term is an oscillating correction 

u (Aq) - u {eq) 
{Ao - eqf 

a^'^ contains the oscillating term coming from the boundaries ±q: 

(635) uT - -^^^l-^tr ([v(-0)(,) , V^-'^-q)] era} . 

Finally, a™ corresponds to the first non-oscillating corrections issued from the endpoints ±q: 

(6.36) < = ^ Z ""^^^^ [v^^'Heq) - [v(^'">(6^)]'} era + u'{eq) tr {[ V(^'l^]'(e?) - IV^^'^eq) [v^''''^]{eq)] cr^ 

e=± 

Here, we precise that u' = dxu, u" = d\u and [y(^'«)]' = dxV^^'"\ 

It now remains to insert the expUcit expressions for y^*-*) as well as so as to obtain the expressions for the 
coefiicients ak,k- 1,2. 

We get that, independently of the time-hke or space-Uke regime, 

"2'" .r . . . ^^ v(q)v(-q) (S- S+ 

(6.37) - i[u{q)- u {-q)\ 



' ' u'{q)u'{-q){2qx)^\S+ S- 

As for ai, we have 

^ox '^1 1 f - / ^ [u{-q)-u{AQ)\SQ . [u{q) - u{Aq)] Sq 
(6.38) — = — j\iv{-q)-— — — -^-^-iv{q) 



xi 2y[^h{AQ)x^ I u'{-q){AQ + q)'^ S- u' {q){AQ- qf S+ 

in the time-like regime, and 

^r,, ai 1 f . . ^ [" ('^) ~ " (~^)] ^- ■ , . {u{M)-u iq)] S+ 

(6.39) -r ^ rVv(-^)-— — — -5-t;--«v(^) 



xi 2 ^/nh (Aq) xi I u'(-q) (Ao + q)^ 'So u'(q) (Ao - qf Sq 

in the space-Uke one. We remind that S± and Sq have been defined in (2.7) and (2.8) 



• The constant term 



The x-derivative cannot fix the constant in x part of the leading asymptotics. We use the /Iq -derivative identity so 
as to fix the Aq dependent part of this constant. Then, in the space-like regime one obtains the io-independent part 
of the constant term by sending Aq ^ co (the asymptotic expansion is uniform in Aq lying uniformly away to the 
right from q). In such a limit, the determinant can be related, up to O (x~°°) corrections, to the generahzed sine 
kernel determiant studied in [32]. In this way, we are able to fix the constant in this regime. In the time-like regime, 
in order to fully fix the constant, one has also to compute the ^-derivative of the determinant asymptotically. 



34 



RIEMANN-HILBERT APPROACH TO THE TIME-DEPENDENT GENERALIZED SINE KERNEL. 



We already know from the above analysis that log det [I + V]- xa-i + log xao + C[v,u,g\ + o(l). Using (3.9), 
we get 

(6.40) [y, ^] = ^ ^ idA,u) (z) |] ^ , = -5^ {v" iq) log |m' (^)| + (-q) log m' (-^)} . 

y(0) 

The absolute value has been chosen so as to treat the space-like and time-like regimes simultaneously. 

In the space-like regime, the asymptotics are uniform in Aq, as long as Aq remains uniformly away from q. 
Hence, one can set Ao = oo in the asymptotics so as to fix the constant term. When Aq = +oo, the function u has 
no saddle-point, a straightforward computation shows that V iA,/u) - Vgsk iA,iJ.) + O (x~°°), with 



VasK iA,M) = -{l- e^'-^^f {l - c^'^^}^ ■ '"'^^^^^^jf^f '^^^ with 7(A) = e (A) (e-^'-(^> - l)' 
Moreover, the big O symbol is uniform on [ ; ^ ]. This means that 



2 



(6.41) det [I + V]= det [/ + Vgsk] • (1 + O (x"")) . 

[-<i'i] 

This last identity stems from the fact that the resolvent of a generalized sine kernel is polynomially bounded in x, 
and this uniformly on [-q ; q], cf [32]. Using the x +oo asymptotic behavior of det|^_^.^j [/ + Vgsk] obtained 
in [32], we get that 



(6.42) C [y, u, g] = -y" (q) log [2q (u'iq) + iO*)] - y" (-q) log [2qu' (-q)] + log G (1, y (q)) G (1, y (-q)) + Ci [v] 

-1 



q q 

J dAg'{A)v{A) - J d^y(^)log'(e-2"''(^)- l) . 



The functional Ci has been defined in (2.4) and we agree upon the shorthand notation G (1, z) = G (1 + z) G (1 - z) 
for the product of two Barnes functions. Note that the iO^ regularization only matters in the time-like regime where 
u' (q) < 0. Of course, for the moment we have only proven the value of the constant term in the space-like regime. 
To see that the constant term is indeed given by C [v, u,g] (6.42) in the time-like regime as well, we apply the 
so-called q derivative method [32]. Namely, starting from the identity 

(6.43) dg log det [I + V]^R {q, q) + R {-q, -q) , 

one replaces the resolvent R by its leading in x part corresponding to sending IT = /2 in the reconstruction formula 
for I (/I) ^ in terms of x- The leading resolvent around ±17 is then expressed in terms of CHF with the use 
of identities (A.6). Then, following word for word the steps described in [32] one obtains that, in the time-like 
regime, dqC [y, u, g] is indeed given by the partial ^-derivative of (6.42). This fixes the Aq and q dependent part of 
the constant term in this regime. As the remaining Aq and q independent part has to be the same in both regimes, 
the constant term is fully fixed. 

The form of the asymptotic expansion given in theorem 2. 1 follows once upon applying the identity 
(6.44) 

- Jlog'(e-2"''(-*>-l)v(i)di ^, , ^ 

e -0 G(l,y(^))G(l,y(-^)) = e'5('^«)-'^-^))(2;r)''(-«)-''^^^G2(l +y(^))G2(l-y(-^)) . 



This identity is a direct consequence of (A.8). 
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7 Natte Series for the determinant 

In this section, we derive a new series representation, that we call the Natte series, for det [I + V]. Just as a 
Fredholm series is well adapted for computing the determinant of the operator I + V perturbatively when the 
kernel V is small, the Natte series is built in such a way that it is immediately fit for an asymptotic analysis of the 
determinant. The form and existence of the series is closely related to the fact that the asymptotic behavior of this 
determinant can be obtained by an application of the Deift-Zhou steepest-descent method. 

Let / + A be the jump matrix for 11. Then, according to sections 3 and 5, A has an asymptotic expansion that is 
valid uniformly on the contour Sn: 



(7-1) -here 



n>0 



A("^(z; x) = ^|x- -3-70- for z e 80x^,6 

(z-^o) 

(n + l)!(z + ^ri ^' 



and everywhere else A*^"^ (z; x) - 0. In other words A (z; x) is a O (x °°) everywhere else on the contour. Moreover, 
one can convmce oneself that this O (x-°°) holds in the n L'^C'T) sense, for any curve that is lying sufficiently 
close to '^E- Finally, we remind that cr = cr+ in the space-like regime and cr = cr" in the time-like regime. 



7.1 The leading Natte series 

We start the derivation of the Natte series by providing a convenient integral representation for log det [/ + V]. 

Lemma 7.1. Let V be the kernel defined in (l.I) and n(/l) = 11 (/I; x) be the unique solution to the associated 
RHP. Then, the logarithm of the Fredholm determinant admits the below representation 

{12) log det [/ + V] [y, u, g] - log det [/ + Vf^ [v, u, g] + log det [/ + Vt'''°^ [v, u, g] 

where 

(7.3) log det [/ + Vf^ [v, u, g\ = ix J u'(A) v {X) &X - {y\q) + v\-q)) log x -I- C [y, u, g\ 

-1 

and C [y, u, g\ has been defined in (6.42). Also 
(7.4) 

logdet[/ + y](^"'')[y,M,g] = J dx'^^G(^)tr[5,n(i;x')cr(i)n-i(i;x')] + ^ J ^^~^~^^r^~^ 

+00 y y dD 

The contour y is as defined in Fig. 14 and dD - -dDq^s U —dO-q^s U -dO^^s- Thefiinction G has been defined 
in (6.30), A^^-* in (7.1) and we remind that cr (A) — o-^ly,o)iA) + cr"'"ly(s)uy(i)(/l). 

The convergence of this integral representation is part of the conclusion of the lemma. 



Proof — 
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The formula for the x-derivative of the determinant (3.8) is the starting point of the proof. By re-ordering the 
terms we get, exactly as in the proof of theorem 2.1, 

(7.5) d, log det [/ + V] [v, u, g] = log det [/ + Vf\v, u,g]+R, 

in which 



(7.6) n = -i— 

or] 



77=0+ 



Here the matrix A'^^^ appears in (7.6) as its contribution has already been taken into account in dx log det [/ + V^^^K 
It had thus to be substracted. 

Now, performing exactly the same steps as in the proof of theorem 2.1, we recast the integral in such a way 
that the jy-derivative can be moved inside of the integration symbol. Note that the operation of squeezing the 
contour T{'^e) in (7.6) to y does not affect the term coming from A*^^^ as it is holomorphic outside of dD. Once 
that the //-derivative has been computed, the result follows by an x-integration. This integration is licit as, due to 
the presence of A^^\ the 7/-differentiated integrand behaves as o(logx/x^) , for x ^ +oo, in what concerns the 
non-oscillating contributions and as e'^^'x"'^ for the oscillating terms. Here, v and w are constants such that v e R 
and % {w) > 0. The oscillating contributions are thus also integrable, at least in the Riemann-sense. ■ 

We are now in position to derive the logarithmic Natte series representation for log det [/ -i- V]. 

Theorem 7.1. There exists a sequence of functionals !Fjv [v, u, e^] (x), such that 

il.l) log det [/ + V] [y, u, g\ = log det [/ + Vf^ [y, u,g\ + Yj • 

iV>l 

There exists a positive N-independent constant m (x) such that [v, u, e^] (x)| < [m (x)]^. m (x) is such that 
m (x) = O (x""') where, for 5 > hut small enough 

(7.8) w = ^ min^l/2, 1 - w - 2max |2l [y(e^)]|j and w = 2max| sup |2l [v - v(e^)]| i 
The junctionals Tn \y, u, e*] (x) admit the integral representation 

(7.9) r;v[v,«,e^](x) = |; Z f^^'cf r Tt^^ •M'^'f'^'?=i'^''<^^^^^^^^ 

^€{±1,0} ' mi) 

in terms of the auxiliary functionals 

-G {A) Dn^,(a, {Zjfj^^,x; {6,}^^i) [y, u\ 



(7.10) HN,r(A, {Zj]%,x; {ejfj^^ [y, m] - 



r-\ \ ^ I \ 

{A-zi) {A-Zr+i)Y\[zp-Zp+i) n \Zp-Zp+l) 

p=l p=r+l 



The first summation in (7.9) runs through all possible choices of the variables ek e {+1, 0) subject to the constraint 
2 6jt = 0. Then, one sums over integrals running over the upper/lower part y^^^ of the contour y and also over the 

associated inslotted contour jlijy I 
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For N >2, the Junctionals Di^ are defined as the fimctionals ofv and u that appear in the expansion of 

N 

(7.11) tr[A(z,)...A(zi)crU)V(z,+i)...V(ziv)] = J] DN,rU,{Zj}^=l,x■,{ej}^^^[v,u] e>'--''''' , 

Set=() 
et€{±l,0) 

into different powers ofe^^^^\ For N - I one has 

(7.12) Di(^Zi,x)=tr[(A(zi)-x-iA(i\zi))(r3] V)(^) . 

Above, V is the adjugate matrix to A; V = Comat [A]', so that /a + V corresponds to the jump matrix for 
Finally, -y^^^ denotes that part of the curve y which lies above/below ofLji- Let Pi, P2 stand for the two intersection 
points between y and Sn, cfFig. 14. Then, 1.^^ is the contour equal eveywhere to Sn except in a small vicinity of 

the points Pk, where it avoids these points by below/above. Then \ is realized as the Carthesian product 
of two inslotted contours of length r and N - r: {sjj''-'-}^'^'^'' = {^n"''}^'^^ X {sj/"''}^^ '^^ ■ 

Starting from the definition 6.1 of the matrices A and V on inslotted contours one defines the matrices A 

t/ L 

and V on E^'^ [zk\ as the analytic continuations of A from Sn {zk\- 



Proof — 

The functional Tn {v,u,g\ will be constructed by merging the integral representation (7.4) with the Neumann 
series for 11 (6.1) and 11"^ (6.2). These series converge uniformly in A (and every finite order i-derivative) on 
every open set O such that d(0, Sn) > 0. However, in (7.4), one integrates tr \dx^ (/I) o" {X) (/{)| with a weigth 
along y, where the contour 7 is depicted in Fig. 14. The latter contour intersects Sn- Hence, it contains points 
that are not uniformly away from £n- However, we have already argued that tr [5^n± (/l)cr(/l)n~^(/l)] can be 
analytically continued to a small neighborhood of Sn located to the right/left of Sn. Such an analytic continuation 
can be also performed on the level of the Neumann series for n*^ 

In order to have a Neumann series representation for 11 or H"^ that is uniformly convergent in /I e y^l^ \ {Pi , P2}, 
we use the local analyticity of the jump matrices for IT so as to deform the original jump contour Xn appearing in 
the RHP for H and TT^ (6.1)-(6.2) into the contour eJ/-"-; 



n <-«) = '^-Z J^s^lri) K«} ['^1 « A fe) and n- w) = 2 J^sf^' W K-) I'^l W 

n 



According to these formulae 11, H ^ are holomorphic on some small vicinity of y^l^. However, we do stress that 
these analytic continuation from above and below P^ differ at P^. Hence, for A e y^^^ \ {Pi , P2}, we get 

(7.13) tr [^;^n {X) cr {X) TT^iX)] = J] /iv U, x) 

N>\ 
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with 



(7.14) fyiA. 



N-l ^ 



-dzdy 



infiA-zf(A-y) 



tr 



C^m JUi] (z) A (z, X) o- (A) V (y, x) {'C^,,, j"^ '[A] (y) 



+ 



I 



-dz 



2ot(/1-z)^ 



■tr 



{c^mf[A]a)A(z,;c)cr(^) 



Above, we have insisted on the dependence on x of the matrices A and V. The representation (7.13) allows us to 
define the functional Tn [v^ u, e^] : 



(7.15) 

and, for > 2, 
(7.16) 



r r\ 1 f dz trlA(')(z;x')(r3l 

+00 y y do 

X 

Tn [v, u, e«] = Jdx'^fy (A, x') G (A, x') ^ . 



+00 y 

We remind that G is given by (6.30) and above, we have expUcitly insisted on its x-dependence. 

In the following we justify that (7.15)-(7.16) are well defined and that one can exchange the integrals over y 
and [ x ; +oo [ in (7.4) with the summation (7.13). Then, we provide explicit bounds for T^n [v, u, e^] and finally 
outline the steps leading to the derivation of the representation (7.9) for "Fn- 



Exchange of symbols 

Building on the identities: 

|tr (AB)\ < maxj,k (\Bjk\) J] \Ajk\ and W^'^^hhm ' Z WMIlH'^ WBekWiH'^) < 4 ||A||i2(^) \\B\\l2^^^ , 

j,k j,k j,k,f 



and after some algebra one obtains that 



k=23 



T=T4 '•=1 



V ( ^ \N-l-r 



A_ 
Ik 



In the intermediate calculation we have used ||cr(/i)||^c»(^) = 1. By using the estimates (6.4)-(6.5), one gets that for 
T =T or J, 

{'^^nF'^i^(z^)- ^ W^W^-i-^j' ¥k ^ (^n)r ) IIA|l..(zj,) 



and 
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Also, one has 
(7.17) 



<2||A|| 



The estimates and asymptotic expansions of A ensure that there exists an x-independent constant C2 such that 

Hence, 
(7.19) 

2N 

- max id "It/' X:, I! I ^1 In'' + r<Y.Ti)\ with 



(7.18) max {||A||i..(sr) + mii^i^rj < 



c(Sn) 



where c (Xn) = max c (Sn ) 

T€(T,il ^ 



ll/ivl 



< 



^'"(^)-;rc(5:n);€(uV 
k=2,3 



max (d-*=(/,S[[)) [^j [n'^ + c (Ln)) with w = | min(l/2, 1 - w - 2max± |!Ry(±^)|) 



It follows that for x large enough and for N > Nq (with A^qw > 1) {A, x') T,p=No ^ /p bounded 
on y X ] ;ic ; +cx) [ by an integrable function. The terms corresponding to p = 1, . . . , A'^o - 1 are also integrable as 
will be shown below. Hence, by the dominated convergence theorem one can exchange the summation and the 
integration symbols leading to (7.7). It now remains to provide sharper estimates for each summand. 

Sharper estimates for Tn [v, u, e^] 

It follows from proposition C.l apphed to the jump contour that for A e 7^^^ one has the representation 



[N/2] h [NI2]-h 

XIn{A)^An{A) + YjTj Z 

h=Q p=Q m=h-[N 12] 



t {q\ x) 
e {—q\ x) 



m-rjp 



(:{Aq;x) 
e {-q; x) 



rjh 



X [e(vox)]¥Hj!f''^(^)[e(v,;x)]- 

ee{±l,0| 



2 



where t{q,x) = Qi^"^i) x-^''^'i\ t{-q,x) = e^>'^-<i) x^^^-i'^ and t{Ao,x) = e'^''(^)x-i. Also, we remind that 1/ = 1 in 
the space-like regmie and 17 = -1 in the time-Uke and we agree upon vq - Aq and v+q = ±q. The matrix A^iA) 
contains exponentially small corrections in x and the remaining part represents the algebraically small ones. 



This representation ensures that for Aey 

. , [N/2] + l h [NI2]+\-b 

UN {A, X) 



(7.20) 



/n (A, x) = 



{A,x) , J_ ■y' ■y' ■y' 



to Um=b^l2]-M-'l^''^l V{-q,x) 



-N 



(A,x) 



The functions c^'* ''^ (A, x) and (A, x) can be expressed as traces involving appropriate combinations of the 
matrices and We have included all the exponentially small corrections stemming from the Afs, 

j - l,...,N into On (A, x). 

It follows from the properties of Aj(A) and n^f^'^'^^ (A), j - I,. . .,N that these functions are smooth in A e y 

and X. Moreover, by using the estimates for the L°° norms of the aforementioned matrices A^ and H^^^'^^ (C.4), 
after some algebra one shows that, for x-large enough, given any k e N there exists an A^-independent constant 
C > such that 



(7.21) 



|aiv(/l,x)| < and 



■N 



(A,x) 



< C^x^^ , uniformly in A e y^'^ 



These estimates remain unchanged when considering first order partial derivatives in respect to x of these func- 
tions. Hence for all integers m, b, p of interest the function 



(7.22) 



-N 



e(^,y) Y''' (m,b,p) 



iA,y)G(A;y) 
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is Riemaim-integrable on 7 X ] x ; +00 [. Suppose that m,b or p is non-zero. Then, for > 2 an integration by 
parts leads to the estimate 

+00 



(7.23) 



[cf 



with w being defined as in (7.19). When m - b - p - one simply deals with a non-oscillating integral. In that 
case, 

[Cf 



(7.24) 



+ 00 



There are two cases of interest to consider. If w = 3/8, then since w = 0(6), taking 6 sufficiently small we get 
that Nw < N (I -w) - I. It remains to treat the case when, for all 5 > small enough w < 3/8. In other words, 
1 - w - 2max± |!!ly(±^)| < 1/2. Therefore 



(7.25) 



1 w I I 
< max hRy(±<7) . 



Thus, taking 6 small enough, so that w < 1/10 one gets max+ |5ly(±^)| > 1/5. Hence, for > 2 



(7.26) 



-1-00 



[C] 



N 



y-NW 



\G{A;x)\ . 



Thus, once upon the integration over /I e 7^ U 7^ = 7, we get that there exists a constant m{x) -O (x~^) such that 
|!F/v [y, u, g\{x)\ < [m (x)f for N > 2. The fact that \ fi [v, u, g] (x)| < m (x) from a direct calculation based on the 
representation (7.15) and the first few terms of the asymptotic expansion of the matrix A. 



Justification of (7.9) 

We conclude this proof by explaining how one can obtain a slightly more convenient representation for each 
individual functional 'Fn [>'. u,g] (x). Starting from the Neumann series representations (6.8) and (6.9), it follows 
that for /I e 7T/J- 

JL r aN„ 



(7.27) tr[5^n(i)cr(^)n-i(^)] = -2 J] f 



'z tr [A (-,) . . . A (- 1 ) (r (.1) V (-,+ 1 ) ... V (-,v )] 



t/ilfr.'V) (^-^l) (A - Zr+l) 11 (Zp - Zp+l) n \Zp-Zp+l) 

I p=l ^ ' p=r+l ^ ' 

Above, {sJ/Y'^^^ is Cartesian product of two inslotted contours {sj/"'-}^'^^ x {s]/-'-}^^"'^.^ To obtain (7.27) it is 
enough to multiply out the two series for 5^11 (/I), 11"^ (/I) and then take the trace. One can convice oneself that 
the matrices A and V are such that 

g(z)g'3 g(;)g'3 g(z)g-3 8(2)0-3 

(7.28) A (z) = e~ A|g=o (z) e ~ and V (z) = e~ V|g=o iz) e ~ . 
This means that there exists functionals D^j^, {Zjl^j, x; {£7)^1) [v, u\ such that 

(7.29) tr [A (zr) . . . A (zi) cr (^) V (z.+i) ... V (zn)] = ^^.^(^' ^^J^U ' {^^^=i) t^' { Z i^p) I • 

Set=0 p=l 
e*6{±l,0) 

Above, the sum is taken over all possible choices of A'^ integers et € {±1,0), such that T,k=i - 0. We replace 
the trace in (7.27) by (7.29) and then insert the result into the integral representation for logdet [/ + Vf^^. One 
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can exchange the 7 and y-integrals with the summation over in virtue of the previous discussion. One can 
also pull-out the finite sum over e {±1,0} out of the integrals. Indeed, given any choice of {ejt), the function 
(A, x') G (A, x') Dj^r is Riemann-integrable along y x ] x ; +00 [ . This stems from the fact that D/v,. is bounded 
in A, and as follows from the previous discussion, is at least Riemann-integrable in x as an oscillatory integral. 
Moreover, it is clear that by harping on the steps that allow one to prove the expansion for Yi^ given in proposition 
6. 1 , one can just as well prove a similar type of expansion for Tn [v, m, e^] . ■ 

Lemma 7.2, Under the assumptions of section 2.1, the Fredholm determinant of I + V, with V being given by 
(1.1), admits the below, absolutely convergent for x-large enough, Natte series representation: 

(7.30) det [I + V] [y, u, g] = det [/ + Vf^ [v, u, g] U + J] J] Z ^« ' ^^'^'^^w) ^tej,,/"^"^] ] ■ 

^ n>l k€'Kn{e,]€S(k) ' 



The term of this series is a O ([m (x)]""), with m{x) - O {x^^) and w being given as in (7.8). The functional 
H„ appearing above is a linear functional in respect to thefuntion Wt^j^^^ Q<^tgizt) Q^f]^ n-variables Zt- It produces 
a weighted integration of this function over curves lying in some small neighborhood of the real axis: 



(7.31) ^„(x,{.,,{..w,j[v,.ne-^-]^2: 2: n S^4'-t / ^ 

.=i..=m J [^p^^ 

^ fl FT • n ^'^-r^Aa, {Zdj]%,, xa; {e,,j}%,) [v, u\ ■ f] e^'^^^') . 
In (7.30), the sum is carried out over all the possible choices ofn-uples of integers k - (ki, . . . ,k„) belonging to 

(7.32) 'K„^l^k^(ku...,kn)eN" : |] ^'^^ = «| 
Each such n-uple of integers defines a set of triplets 

(7.33) Jik] = [{s,p,j) , sell; nj, pell; k,^, jell; s]i] 
and a set of doublets 

(7.34) Dik] = [{s,p) , sell;nj, pell;k,]\] . 

A triplet {s, p, j) belonging to J[k] is denoted by t and a doublet {s, p) belonging to Dj^^j is denoted by d - (s, p). 
The notation d\ stands for the first coordinate of d, ie if d = is,p), then s - d\. Once that a choice of k is made, 
one sums over all the possible elements of 

(7.35) eW^jteW,,, : ^£{±1,0} and YA^d,j^O for all e . 



The sums and integrations in (7.31) are also ordered by the sets of triplets J[k] and doublets One first starts 
to sum up over r^, where d runs through D[k]- There are #D[k] such sums in total, corresponding to d running 
through the set D[k]- Finally, for each d e D[k], there is one integral over the corresponding x^, one over the 

corresponding Aa and d\ integrals over the subordinate set of z-variables \zd,j\ .' • 



Proof — 
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We define an auxiliary function 
(7.36) A (r) = 2 /"Tn [v, u, g] (x) 

N>1 

which is holomorphic in 7 on the open disc of radius m~^(x), with m(x) = O (x^^) and w given by (7.8), as follows 
from the estimates on the growth of Tn with x. This means that e'^^^^ is holomorphic on the same disc. The radius 
of convergence of its Taylor series around 7 = has its lower bound given by (x). 

The series for the Fredholm determinant is obtained by using the Faa-di-Bruno formula so as to compute the 
derivative of e'*^''^ at 7 = 
(7.37) 



1 ^ 

nl dy" 



Air) 



1 d" 
nl dy 



Air) 



y=0 



[4m (x)]" d" 

< ; ■ — exp 

nl dy" ^ 



1 



2-7 



y=0 



I.sks=n 



Where we used that \ fM [y,u,g]\ ^ itn{x))^, with m{x) = 0(x The last estimates allow one to see explicitly 
that the Taylor series at 7 = for e'*^'*^ has a radius of convergence that scales as (x). It is in particular 
convergent at 7 = 1 leading to 



(7.38) det[/ + y][y,M,g] = det[/ + y](o^[y,M,g] 1+2 ^ FK^^O ' ]"[ l^^' t^' ^^^1 

1^ n>l 'Lsks=n j=l j=l 

In this language of doublets and triplets, the expression for the product in (7.38) reads 

(7.39, fi(r,i.u.su4' - 1 z zn i^^t / 



7=1 

eij€{±l,0} 



X Y\ ^dur4(^dAzdj]%i,Xd;{€dJ''jli)[v,u] exp] ^t8izt)\ 



d€Dfk) 



The result follows. 



The expression (7.31) for the functional involved in the Natte series is more explicit then as it was given in 
theorem 2.2. 



7.2 Proof of theorem 2.2 



The first part of theorem 2.2, ie the very form of the expansion (2.9) is a consequence of lemma 7.2. The latter 
provides moreover a more exphcit form for the functional 'Hn [v, e^,u]. 

The well ordered asymptotic expansion in x for each functional [v, e^, m], as given in (2.12), is a di- 
rect consequence of the existence of a similar representation for f^^ [v, u,g] together the correspondence (7.38) 
between !f]v[v, M,g] and det[/+ V]. Finally, the existence of a representation for 'J^i^[v,u,g] in the spirit of 
(2.12) can be readily obtained by inserting the well-ordered series representation for 11 and (we remind that 
11"^ = 'Comat(n) since det [11] = 1) given in proposition 6.1 into the integral representation for 'Fn, (7.14), 
(7.15), (7.16). ■ 
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7.3 Higher order Natte series 



The higher order Natte series is a generalization of the Natte series derived in the previous sub-sections. It gives a 
direct access to part of the asymptotic expansion without having to compute the effective form of the functionals 
Dyv.r and then carry out the contour integrals. Indeed, even if it is possible in principle to compute explicitly, order- 
by-order the functionals D^/^r and thus H„, this task becomes very quickly monstrously cumbersome. In order to 
get the corrections, it is more desirable to apply the procedure below (or its obvious extension to higher order 
asymptotics) if one wants to access to the higher order correction then those contained in det [/ + V]^^^ [v, u, g]. 

Proposition 7.1. The Fredholm determinant ofl + V admits the below convergent Natte series representation: 
(7.40) det [/ + V] [v, u, g\ = det [/ + Vf ^ [v, u, g\ exp | ^ &x' [x']"'^ ai{x') + [x']'^ [a^''(x') + 4°(x')] 

V+oo J 

+ Z ^.(^,W.k*})h«,n,,,,^,e^'^^^'^] |. 

^ n>l ke'K„ [£,]ea{k) ' 

There //„ is defined as in (7.31), (7.10), but with the minor difference that the functionals in (7.10) should 
be replaced by the fimctionals cis given in (7.42). Also, a\, a^'^, a^^'^ are given by (6.36), (6.37), (6.38), 
(6.39). Note that here we have explicitly insisted on their dependence on the large-parameter xf. The fundamental 
difference between the higher order Natte series and the one discussed previously is that for x large enough and 
for an n-indepenent O; 



HnixAk) ,{et)) 



• 1^1 M ailogx a°''''logx ai ^2 

= min<^0 — , O' 



5 
X2 



X2 



■2+w 



The constant w is as defined in (7.8). 



Proof — One starts by performing the decomposition 

X 

log det [/ + y] [y, u, g\ = log det [/ + Vf^ [y, m, g]+ J dy 



(aiiy) al''{y) + a^^{y)^ 



3 



+ 



+\ogAQi[I+vt^''^\v,u,g\ 



There, log det [/ + u, g\ corresponds to that part of log det [/ + y](^"''^[y, u, g\ (7.4) where all terms that 

give rise to the integral involving ai and a°2^^^° have been substracted. Namely, 

X ^ 

(7.41) log det [/ + y](^"^2) ^ g^^ Jdx' ^ ^ G {A) jtr [d^H {A) a {A) IT^ {A)] - ^:^tr {[n(i>]'(^) cr^] 

+00 y 

{[n(2)]'(^) _ n(o)(^) [n(i)]'(^) - n«(i) [u^'^]\a)] a, - ^tr ([rf3)]'(^) - n«(A) [rfi)]'(^)) ^.3} . 



3 

yl 



We stress that the variable of integration x' corresponds to the large parameter (denoted by x before) that enters in 
the formulation of the RHP for IT and on which IT depends implicitly. 

The expansion of log det [7 + V]^™''^ [y, u,g\ goes along the same lines as before, with the minor difference that 
the functionals Di^r are defined slightly differently. Indeed one has to substract from the Neumann series like 
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expansion for tr [^^II (A) <t (A) IT ^/l)] all the subleading contributions that appear in (7.41). For this purpose we 



define 



(7.42) tr [A izr) . . . A{zi) (r{A)V (z,+i) ... V (zn)] - %r {{z} ,x)^ ^ £>^,, (a, [zjf^ , x; [ejf. ^ ) [v, u] 

e*€{±l,0) 

In order to define R^^^r ({z} , x) we represent the asymptotic expansion of A (z) sUghtly differently then in (7.1). 

< (fP\z) 



2 epg(zp) 



(7.43) A (z) ^ 2 x" 5 (z, X) with 

p>i _ ^^^^^^ ^ 

£=±p\{z-eqy 

There we took 6 small enough and 5 > 6' > 0.\n terms of such matrices one has 



y(e;p-l)/- \ 

M(2^)(z;x) = Z-^ ^.i^^^^^^^^^^^Xz) 



^iV,r(k},Jc) = (5iv,i(5;., 



N 



2x-^tr[M(^>(zi,x)cr3] 
^=1 



+ 



tr [m(1)(z3, x) m(2)(z2, x) m(i>(zi, x) cr^] 



r=l 



2 x-^tT[M^J'\zux)(r3M^P'\z2,x)\ 

p,p'=l 



Finally, defining //„ as in (7.31), (7.10), but with the minor difference that the functionals Df^ in (7.10) should be 
replaced by the functionals Dn given in (7.42). One gets the desired representation. ■ 

A similar Natte series can be obtained for other quantities that are also related with the correlation functions in 
integrable models. 

Proposition 7.2. Let F\ be as defined in (3.3), then the below Fredholm minor admits a Natte series representa- 
tion: 



(7.44) 



■sin^[7rv(^)] Fi {A)E{A) 



det [/ + y] [y, u, g] = det [/ + Vf^ [v, u, g\ 



r&A fdi^ . 

V^wW ^«'(^) " ;c«'(-^)'^- ^ 4^ 2. Z (x,*,te})[y,.,n,,,,^,e J I . 

r^ere summation runs through all the possible choice of integers k\,... km+i belonging to 



(7.45) 



'Kn = \k ^ {ki,. . .,kn+i) : ks en , s ^ 1,. . .,n and fc„+i e 



kn+i +Y^sks = n[. 

s=l ' 



The remaining summation run through sets that are labelled by doublets and triplets belonging to 
J{k] = {(s,p,j),se\[l;n+lj, pell; ks^, jell; s-ds,n+in]li] 
D{k} = [(s,p,j),sell;n'i, pell;kslll} . 

Indeed, then 

(7.46) &{k)^\ {^t}t^Jm ■ e {±1,0} , Yfd,j = with d e D^k) and J] e^t^^i.^i = 1 \ 
y j=i p=i f 
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Finally, the functionals {x, {k} , {e<}) read 



v,u, Y\ e^'^^"^ 



n 



2in 



when k„+i — n and in all other cases, 
H^:\xAk},{et}t^j,,^ 



2in 



r ^7' 



Afe?e that we have used above the notation introduced in lemma 7.2. Also we have set 

2 



p=i 



i- 



N-l 

n [zp - zp-i 



) 



CnHzj}^, X, {ej}^) [y, u] f[ e^''^^^") 
p=i 



W-1 



itr [A (zjv) . . . A (zi) o- ] • ]~[ (zp - Zp-i) ^ 



J] C;v({z,}f , X, ) [y, «] ]~[ M^p) 

p=i 



r/ie 5Mm5^ in the two equations above run over all choices of the variables e^, s - I, . . . ,N with 6, e {±1,0} and 
e = 1 



Proof - 



By using the integral representation forx (3.6), one readily gets that 



(7.47) 



— 4sin2[;ry(i)]Fi (^)£ 
2;r 



(i) 



Also, it is easy to convince oneself that 

(7.48) i [no.]i2 = i \xoo]n + J ^"'(^) ^ ■ 

The claim follows by expanding [IlooJia into a higher order Natte series (where the first few terms of the asymp- 
totics have been taken into account) and then taking the product of this series with the Natte series for the Fredholm 
determinant. The details are left to the reader. ■ 

It follows from the leading asymptotics given in (7.44) that, in the case of the time-like regime, the saddle-point 
Ao does not contribute to the leading order. It can however be checked that it does eventually contribute. Its 
contribution is a O^x'^'^, x'^'^ix-''^'''' + x'^^-i^f). 
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8 Conclusion 



In this paper we have obtained the first few terms in the leading x +00 asymptotics of the Fredholm determinant 
of a class of integrable integral operator that provide a starting point for the analysis of the large-distance/long- 
time asymptotic behavior in integrable models away from their free fermion point. Also, we have derived a new 
series representation for the Fredholm determinant, the so-called Natte series. This series is well adapted for an 
asymptotic analysis of the Fredholm determinant and can thus be thought of as being an analogue of the Mellin- 
Bames integral representation for hypergeometric functions. In two subsequent paper, the Natte series will appear 
as a central tool in computing the large-distance/long-time asymptotic behavior of the correlation functions in the 
non-Unear Schrodinger model away from its free fermion point [34, 37]. As a byproduct of our analysis, we have 
been able to bring a little more order to the structure of the asymptotic expansion of Fredholm determinants of 
operators belonging to the class of the generalized sine kernel. It would be interesting to extend/push forward the 
form of the full asymptotic expansion of the determinant given in theorem 2.2, in particular by providing a closed 
form (ie the explicit values of coefficients/functionals), at least in the case of some particular kernel such as the 
sine kernel. 
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Appendix A Several Properties of CHF 



One can check that for z e R"^ 
(A.l) 



^|l,-;-e''^j-Y|l,-;-e-zj-2/^-e--. 



*F (a, c; z) has an asymptotic expansion at z ^ 00 given by 
(A.2) y(a,c;z)^|;(-ir ^"^"^"-;^^^" z— -hO(z-^-) , -| 

n=0 

It also satisfies to the monodromy properties 
(A.3) 

(A.4) 



r^(a) 

^(a, Uze-^"") = >F(a, l;^)^^'™ - y(l - a, 1; -z). 



THa) 



< arg(z) < y . 



< 0, 
3(z) > 0. 



Tricomi's CHF can be expressed in terms of Humbert's CHF 



(A.5) 



^(a,c;z) = r 



1 -c 
a- c + I ^ 



cD(a,c;z) + r 



c- 1 



z^~''<b(a-c+ 1,2- c;z) 
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There exists a similar formula expression Tricomi's CHF ^ (a, c; z) in terms of Humbert's one 

(A.6) O (a, c; z) = p^^e^^T (a, c; z) + L^Jenia-o^z^ _ ^. ^ 

r (c - a) r (a) 



where e = sgn (3(z)), and 

(A.7) 0(a,c;z) = yi^^fi 

The Barnes' G function satisfies to the reflection property 



ia)n z"_ 
T' 



G(l+z) 



(A.8) G(l -z) = ^^^^^ exp ^ ( ;r^cot[7r;v:] d;c|> , 

which holds for 21 (z) < 1 in the usual sense (and also everywhere else by analytic continuation). 

Appendix B Proof of the asymptotic expansion for n 

B.l Two lemmas 

We first need a technical lemma 

Lemma B.l. Let the matrices Aj take the form 

where the entries aj, bj, cj, dj do not depend on tj. Then 



0-3 



(B.2) A„ . . . A, = X ( A!^ 'd'" ] • ( ^^^^l 

a=0 y(a)gs„.jv ^ I \ n na-l/ 

Above, the sum runs through all choices ofla-uples of integers jf^^ - . ., 72a) with /"^ belonging to 

(B.3) Sa;N = [Uu iia) e [N*]''' : 1 < ;i < • • • < ;2a < ^) • 

The entries of each matrix appearing in the sum are linear polynomials the entries of the matrices A^. 

This lemma allow us to trace back all the dependence on the fractional power of xin the products A {zn) • • • A (zi) 
of the non-trivial parts of the jump matrices for 11. 

Proof — 

The result clearly holds iorN=\ as then, the only possibility is to take a-0. 
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We prove the induction hypothesis for the 1 1 entry. It goes similarly for all the others. By applying the induction 
hypothesis to Ajv . . . Ai and then expUcitly multiplying out with A^v+i, we get that 

[f] 

(B.4) [A;v+i...Ai]ii = y y I ' ' ' \ (ajv+iA .(a) + ^-^bN^iCiaX 



The result follows as the above sums can be seen organized in respect to the partition 

[N+lll] .[/V/2] ^ .[Nil] . 

(B.5) U Sa-^^l^\ U^«;^fUi [j[^^h<---<j2a<N-\, j2a+l^N , j2a^2^N+\]\ 
a=\ ^ a=l ' ^ a=l ' 

r [Nil] . 
U U [^^jl<---<ha-l<N-\, j2^^N+\]\. 
^ a=l ' 

■ 

Lemma B.2. Let 'J^pj (zi, . . . , zn) be a holomorphic function on D = i)v,Si X • • • x Dv.s^, where < 6n < ■ ■ ■ < Si 
and V e C. Let dD = dD^^Si x • • • X dOy^s^f be the skeleton ofD and np a set of positive integers. Then, for A lying 

outside ofOvSv ^'^^ 
(B.6) 

r d^z rN(zu...,ZN) fi 1 _ y y 1 [ 1 g^' ^^^1 

k=\ 

where we agree upon 

N N 

(B.7) rp^Yjf^f + ^-P- no = 0. 

e=p e=p+i 



The proof is a straightforward induction. Note that the total highest possible order of derivatives of Tn that 
is produced by the above contour integral is «^ + - 1. It corresponds to no-derivation in respect to the 
variables zi,---,zn and a derivative in respect to zi of order n^ + N -\. All other choices of the integers {kg] 
lead to a total order of the partial derivatives that is strictly smaller. 



B.2 Proof of proposition 6.1 

We are now in position to prove proposition 6. 1 . 

Here, we will only discuss the case of the time-Uke regime. The proof in the case of the space-Uke regime goes 
very similarly, so that we omit it here. 
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We have already established that, for ;c-large enough, n (A) is given in terms of a uniformly convergent Neu- 
mann series (6.8): 

A(z;v)...A(zi) 



(B.8) n (A) = i2 + Y. ^ = ^'' { 4- 

N>i ^ ± (2?7r 



m {A - zi) n [Zp - Zp+\ ) 

p=i 



Above each N-foId integral runs across the inslotted contour as defined in Fig. 1 1 and the equality holds for A 
uniformly away from the boundary En- We remind that in this Neumann series the matrices A (zk) are subordinate 
to the jump matrix I2 + A (A) for Il{A) solving the Il-type RHP associated with the jump contour Xn [Zk], cf 
definition 6. 1 . 

To prove the claim of proposition 6.1, we build on (B.8) so as to obtain a more precise form of the asymptotic 
expansion of n^v (A). 

Recall that each contour En [Zk] entering in the definition of the inlsotted contour l!^^ can be divided into its 
exterior part T [zk] and three circles dDq^Sk (^^-qA ^2)/io,<5yf There < 6n < ■ ■ ■ < di and 61 is small enough, 
in particular it is such that di < \Aq ± q\ /2. However, the very choice of the contour f implies that 

(B.9) m\L^(^u.^) + ll^llL^(rb.]) + \\^h^(^u.^) = ^ (^-°°) • 

Hence, from the point of view of the asymptotic expansion, one can drop all contributions to (A) stemming 
from those parts of the multiple integral in (B.8), where at least one variable is integrated along F. Indeed, due to 
the estimates (B.9), such an integration would only produce O (x~°°) terms. 

The matrix An(A) appearing in (6.11) contains exactly these contributions, and hence A^iA) - 0(x~°°) uni- 
formly away from Sn. 

It thus now remains to focus on the efi'ect of the integration on the boundary of the three disks centered at ±q 
and Aq. In other words. 



^_d^z A(z;v)...A(zi) 



(B.IO) (A) = An (A) + x" 2^ &) -— ^ • — where zq = A . 

k=\ 

The above sum corresponds to summing up over all the possible choices of the integration contour dDy^^^Sk ^'^^ 
each variables Zk- More precisely, one sums over all the A^-dimensional vectors e belonging to 

(B.ll) £w-|e = (ei,...,6A,) : ek e {±1,0}} . 

We also agree upon the shorthand notation v+ = q,v- = -q and vq = Aq. Finally, the integration contour dDe in 
each summand corresponds to the Cartesian product of N-circles dD^ - dDv^^,Si x • • • x ^Dy^^^s^ of decreasing 
radii < 6j\f < ■ ■ ■ < 61, with 61 small enough. 

We stress that there exists natural constraints on the possible choices of the e^. Indeed, if zj and z^+i both belong 
to a sufficiently small neighborhood of Aq, then A.(zjj ^(^Zj+i^ - 0. Hence, choices of A^-dimensional vectors e 
having two neighboring coordinates (ej and e^+i for some f) equal to zero do not contribute to the sum in (B.IO). 



The asymptotic expansions of A (z) on each of the three disks all take the generic form: 

e)]? •AW(z)-[e(z;e 

n>0 ^ 



... [e(z;e)]^ •AW(z)-[e(z;e)]-^ -f 1 • . ,n \ .n ^ m 

(B.12) Aiz)=^2_, ~. „(2-|g|)+i uniformly in z e £)v„2,5 \ £)v„<5' ee{±l,0} 
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The radii are such that 6 > 6' > and 6 is taken sufficiently small, but are arbitrary otherwise. (B.12) is to be 
understood in the sense of an asymptote expansion, ie up to a truncation to any given order in x. The detailed 
expression for the matrices A*^"^ (z) and e (z; e) differ on each of the disks (ie for e = ±1 or 0). However, e (z; e) 
are holomorphic on any sufficently small neighborhood of ±q or Aq. Also, the matrix A'^"^(z) does not depend on 
X anymore. The function e (z; e) contains a fractional power of x and also an oscillating term: 



(B.13) 



efce) = 



Jxu(-q) Jlviz) 



aixu(Ao) 



X 2 



for e = 1 
for e = -l 
for e = 



We are now in position to establish the asymptotic expansion of the second term in (B.IO). 
Expanding each matrix A(z,j) into its asymptotic series (B.12), using that the latter is uniform on the compact 



domain of integration we obtain the asymptotic expansion of li^: 
fB 14) n a) - V - V V f-^ [e^]? A("~)(zjv) [ejv-i/e^]T. . [ci/e;]^ AC'^zi) [ei]"^ 



r>0 e^SN „e/4> aO 

There appears a summation over A^-dimensional integer valued vectors n belonging to 



N ^ I V' 

U - ^l) • n {Zk-l -Zk) Y{ [Zp - Ve ) 
k=2 p=l ^ ' 



(B.15) - {r = . . . , «w) : «yt e N , ij^i'-,. ■ , - 

Note that in order to lighten the notations slightly, we have set e^t = e {zk\ ^k)- Also, just as in (B.12), we did not 
make the remainder expUcit. 



Lemma B.l ensures 
(B.16) 

[e^]^A(»^)(ziv) 



(n) (ri) 

the existence of holomorphic functions A^^^, {[Zk]), D\^^y {{Zk}) of zi, ... ,Zn such that 





2 











[f] 



A^''^)(zi)[ei]-^= J] J] 



' ^,i\z^ e^B2({z,})) 



^-^C'^faAZk}) D^^iiZk}) 







0-3 









Due to the form taken by the matrices A^"^(z), not all configurations of the 2a-uples /"^ appear in (B.16). 
Indeed, when Zk £ dO\^^Sk - 0) the matrix A^"^(zjt) is proportional to cr~ {cf (5.14)). It appears in (B.14) 

with a pre-f actor e^^. Therefore, for Zk £ dDAo^St the only non-vanishing terms in the sum over f'^^ € Sa;„ are 
those corresponding to choices of 2a-uples f "^ = (ji, . . . , j2a) such that jp - k for some p. In other words, each 
time an integration variable belongs to dD^^^s^ for some k, the associated oscillating exponent t~^(Ao; 0) is always 
present. Moreover, all matrix entries in the expansion (B.16) that appear (after taking the matrix products) in front 

of tiie monomials (e vanish whenever a function e^^, = t[zj/,ej^,) corresponding to 

Zjp e dDAo,dj^ would appear in the numerator. More precisely, if there exists a p such that jp - k then 

• for e 2N + 1 (ie eyjf+i = for some 0, one has = D^"j) = ; 

• for p e 2N (ie ej^, = for some £), one has A^j^^ = = . 

Putting together (B.14) and (B.16) leads to the below form of the asymptotic expansion for IIjv 

(B.17) (^) - Z 17 Z Z Z M ■ 
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There Mj(a) stands for the matrix appearing in the expansion (B.16) and Ie„ is a functional depending on the 
choices of the entries of the A'^-dimensional vectors e and n. It acts on holomorphic functions (or matrices in the 
sense of entry wise action) - BDy^^ x ■ ■ ■ x dOy^^ according to: 



(B.18) Ie;n[rN]-- 

d!D, 



s=2 s=l 



The functional I„-f can be estimated by computing the residues at v^, . This produces partial derivatives of at 
the points = v^^ . From now on, we focus on the analysis of the action of /„;e on the 1 1 entry of the matrix M^<a) . 
The case of all the other entries can be treated similarly. 

Depending on the choice of the components of the N-dimensional vector e and hence of the evaluation points 
Vtp, after performing the integration induced by (and having computed the eventual derivatives) the ratio 

(^72 • • • ^jia) I ("^ii • • • ^ha-\) present in the 11 entry of (B.17) reduces to: 
for some -a < m < a ; 



tP{q; +) t^~P{-q; -) / t{q; +) 
or '— ■ — • — — ■ — - I for some \ < b < m, 0<p<b,-(a-b)<m<a-b; 



e^{-q; -) 
+ 

7*(^o;0) \e(-^;-) 

Hence, we get that there exists two sets of constant c and c^'^'^\ 

h Jt 

(B.19) 

V V 1 W^^ 1 1 V / r (m)^v V V e^(^; +) e^"P(-^; -) / <q\ +) T 

„6//<'-' fe£jv m=-a\^^ >l b=lp=Om=b-a ^V/tO,"; \ei ;/ J 



(m, b, p) 

(a) 



Each derivative of the factor in respect to Zk, when Zk is in a neighborhood of ±q, produces one power of 
log X. This log X term appears due to a dilferentiation of the exponent x^^^*^^^*\ It thus follows that the coefficients 
c^™| and cpl''^^ are polynomials in logx. In the following, we determine the degree of these polynomials. This 
will allow us to show that 

(B.20) maxdeg(c^^^) = r + A^-2m-(5^,o and maxdeg(cp*'''^j = r + - 2(Z? + m) 

where the sup is taken over all possible choices of n, e, Once that (B.20) is estabUshed we get the claim. 



The degree of c\" 



im) 



a) 



As already argued, when - 0, there necessarily appears e (zkl 0) in the denominator of • As no function 

e (zi; 0),£ k, can appear in the numerator, this implies that in such a situation e (Aq, 0) appears with some strictly 
postive exponent after computing the integrals. Therefore, one ends-up with a term that does not corresponds to 
the coefficient c^'"}. Hence, contribution to the coefficients c^T} can only stem from ffiese choice of N-dimensional 

vectors e whose entries are in {±1}. This means ffiat when focusing on c^^,, all "admissible" choices of ffie 
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A^-dimensional vector e can be parameterized as 

(B.21) € = { ei, . . . ,€[ , 62, . . . , 62, . . . , ep,...,ep) with = (-1)'*"^ e\ e\ e {±1} for some p < N . 

We now compute exphcitly the action of the functional Ie-„ corresponding to some e given by (B.21). For this 
purpose, it is convenient to relabel the integration variables appearing in (B. 18) in a form that is subordinate to 
such a representation of the vector e. Namely, 

(B.22) (yi,...,yN)^{zi,i,---,Zi,{i,Z2,u---,Zs,t,---,Zp,tJj, ie Zs,t = y-g^+f , where 4 = ^^^. 



r=l 



We relabel the entries of the vector n in a similar way, ie ns,t - Then, the functional Ie-„ reads 



(B.23) /,;„ [[Mjia,\^] = ^ 



_fz 
(2ot) 



N 



p e, 

nn 



s=i y f=i J 



Here we agree upon €q -Q and zo,o - The above integral is directly computed by an inductive application of 
lemma B. 2: 



(B.24) ^.«[Mn]=Zfin 



1 5^'' 



[Mj4,,{{z))\\{Zs-U..,-Vey 



In (B.24) one sums over integers ks,t with s - \,...,p and t - \,...,£s where each kg^t is summed from to 

(B.25) Ys^t = ^ risj + (s-t - ^ ks,t with n,,o - . 

j=t j=t+i 



It follows that each block of variables {zs,i, ■ ■ •,Zs,(,) associated to the same 6„ is subject to partial derivatives 
of total order Ya/^\ ks,t- Hence, the maximal total order of all the derivatives that may act on this block of variables 
is = " 

s,j + ^.v - 1- The unique way of realizing this maximal order is through a single derivative of order 
with respect to the variable Zs,i- We stress that in this case, all the other variables of the block are simply 
set equal to v^^ . Very similarly, the maximal total order of all the partial derivatives that may act on a sub-block 
of variables {zsj, ■ ■ ■ ,Zs,e,) associated to the same is r^^^ — + - t- The unique way of realizing this 

maximal order is by a derivative of order with respect to the variable Zs,t- Then Zs,t+i, • • ■^Zs,{, should be set 
equal to v^^ . 

As we have already mentioned, the function e^-j = e {yj^ ; etj depends on x. Hence, its derivative in respect to 
yj^ generates powers of log x. Therefore, the highest degree in log x appearing in Ie-„[Mj(a)\ will be generated by 
a derivative of the highest order possible in respect to the variables jj^, with k = I,..., la. 

Hence, to be able to determine this maximal degree in log x, we first have to order the indices jk according to 
the block to which they belong. For this purpose, we set 

(B.26) :^s^{k : ik e E^. + 1 ; ^.+1 1} and a, = min {A: : A: e . 

Suppose that one deals with a block of variables {zs,\, ■ ■ ■ ,Zs,t,) such that #J{s 4^ 0. Then the highest possible 
power of In x that an integration over the variables of this block can produce will be issued by the action of a 
derivative of the highest order possible on the variable . 7 . Thence, an integration over this block of variables 
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generates a polynomial in In x of degree r™^ _ . Clearly, if - 0, its associated set of variables and functions 
cannot generate, once upon being integrated, any power of log x. 

ai - I a2 - 4 

I I T I I ••• I I I I T ••• I I I I I 




J\ Jl J3 J 4 J 5 ]2t-\ ]2( 



Figure 15. Definition of the sets and of its minimal element a^. One has Jii - {1,2,3), 
^2 = {4, 5), ^3 = 0, . . . , = {2^ - 1, li). The eu deUmit the block of variables of length 4 
associated to e^. 



We now characterize the oscillating term appearing in I^ n [^/n)]- If a given block (zs i, . . . ,z.v/,) corresponds 
to a set having an even number of elements e 2N), then after taking the derivatives and once upon 
evaluating at Zs,t = v^^, the associated ratio of the functions tj cancels out. Indeed, there are as many identical 
factors in the denominator that in the numerator. For instance, when Us is even one has 



(B.27) 



= 1 



However, if a given block {zs,i, ■ ■ ■ ,Zs,e,) corresponds to a set Jl,,. with an odd number of elements (#J?l.v e 
+ 1), then after taking the derivatives, the associated ratio of the functions tj reduces to [e (v^^; e^)]*^"^^ Indeed 



(B.28) a, e 



= e {ve/,€s) and e 2N + 1 



1 



e(ve,;e,) 



Therefore, we obtain that 
(B.29) 



[[m^<„,]^ J = P^C'^log x)- Y] [e (ve;, es)]^-'^"' where deg (pI^'^) = J] 



s=l 
#^je2N+l 



s=l 



Now, in order to obtain the coefficient cl^, we should sum up (B.29) over n e N^' and also over all the 

possible configurations of vectors e parameterized as in (B.21) and such that we eventually generate the power 
(e (^7; +) /e (-^7; -))'". Then, among such configurations, we should look for those that correspond to a polynomial 
Pe n (log x) of highest degree. 

Given a fixed number of flips p in (B.21), one maximizes the degree in (B.29) by choosing the lengths ^ is 
such a way that ^ for any s and such that ja, = + 1. One can do so for all s, but 5=1. Indeed, in 
the latter case one necessarily has jai = ji > 1. Therefore, for such a choice of lengths 4, once upon choosing 
«i f = for f = 0, . . . , Jl - 1 one obtains that this maximal degree of is r + - - (ji - 1). Note that, we have 
used Zi,f«i,f = randYjs^s = N. 



There is also a condition on the number of flips that are necessary to generate the oscillatory factors 

(e (q; +) /e (-17; -))*". Due to the form of the oscillatory factor in (B.29), we get that one sequence (e^, . . . , e^) 
generates at most one factor [e {esq; es)Y, t - ±1. Hence, ifm^O there are at least 2m flips necessary to generate 
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the factors (e (q; +) /e (-q; -))'". If m = 0, then one still has one sequence (ei, . . . , ei) of length €i - N. Therefore, 
one has p > max (2m, 1). Taking the lowest possible value, ie p - max (2m, 1), we get that 

(B.30) max deg (p^P^ - n + - 2m - 5^,0 + O'l - 1) • 



In order to obtain estimates for [H^vln one should still sum up over all the possible configurations of 2a-uples 

(m)-] 
N,ri 



f"K Therefore the highest degree in log x of [n^^]jj is obtained by setting That is, we reproduce (B.20) 



The degree of cp*'^^ 

It follows from the previous discussions that each time an integration over ^'Dx^■,^s occurs in (B.17)-(B.18), there 
appears, once upon integrating, a factor ir^{AQ,0). Hence, the oscillating factor in front of c'^.^j '^^ is necessarily 

generated by these choices of A/^-dimensional vectors e where precisely b entries are equal to zero {ie there are 
exactly b integrations over dD^o^s)- 

Taking into account the fact that, as previously argued, two neighboring entries of the vector e cannot simulta- 
neously vanish, we get that such vectors e can be parameterized as 

p 

(B.31) e- i €i, . .. ,ei , . . . ,eri, . . . ,eri,0,€ri+i, . . . ,€r,,,0, . . . , €p,...,ep ) , where ^(r-N-b. 
We relabel the integration variables y,- appearing in (B.18) in a form subordinate to (B.31): 

(B.32) Cyi,...,yw) = {zi,l,...Zr„A„,Wa,ZTa + l,l ' ^^'t^y{,+t' ^"^yira + l 

where we agree upon 

s-l 

(B.33) 'es^Yj^r + #{k : Tk<s} . 

r=l 

We also relabel the entries of the vector n in a similar way, ie ns,t - +, and nf^ = +i- The action of the 
associated functional Ie „ takes the form 



1=1 f=l 

The integrals over cja are readily computed. We set 



(B.35) ^N-bi{z])^Y\ 



Then, the analysis boils down to the case previously studied: 
(B.36) 
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The integrals runs over the contour dDf^ which corresponds to that part of the initial contour dDe where the 
integrals over the variable (Oa have been suppressed. Therefore 



(B.37) 4^4[M^4j=|;f]f[| 1 ^ 



&N-m ({zD 



|Zj,(=Ve, 



.S=l 



The sum over the integers ks,t runs from to r^,; - 21 j' + is-t - Z ks^t • 

j=t ' j=t+i 

Similarly to the previous analysis, we set 
(B.38) = { A: : ^ e [[^^ + 1 ; + 4 J } and ^ mm{ k : k e ] 

It is then easy to see by using similar arguments to those invoked for c^.^| that 
(B.39) 

'•"[K"]..l = 4^-n Wv..;eJl>-»" where deg(p«)=|: 

#^j€2N+l ms*o 
In order to obtain the maximal degree in log x associated to the oscillating term 



^ ns,k + ^s- [ja, - is) 



e*Uo;0) \e(-^;-) 



present in [IlAflii, we should maximize the degree of the previous polynomial in (B.39) under the constraint that 
the sequence in (B . 3 1 ) ought to change its value at least b + m times (this in order to produce the sought form of 
the oscillatory term with its associated power-law behavior) and that these changes are such that eventually (B.40) 
is generated. 

We should also maximize this degree in respect to all the possible choices of 2fl-uples f"^ of various lengths 
2a and over the allowed vectors n € Nf^ ■ In order to obtain this maximal degree, one should choose a minimal 
number of flips (m + b), choose the lengths and the jk in such a way that ja^ = + 1. Finally, one should also 
take nf'' = for all a. This leads to the conclusion that the maximal degree in\nxisr + N -2{m + b). ■ 



Appendix C Fine bounds on n 



N 



In this appendix we provide bounds for the matrices 11^ 
sition 6.1. 



(m,h,p) 



entering in the decomposition for fl^v given in propo- 



Proposition C.l. Let Xn be a contour appearing in the RHP for 11 and U any open set such that d (f/, Hn) > 0. 
Let Hn be as defined by (6.10) and, agreeing upon rj — \ in the space-like regime and t} — -\ in the time-like 
regime, let 



(C.l) 



' ^ixu{q) ^-2v{z) for 6=1 

e(z;e) = \ e'^^i-'ih^vij^ for e = -1 
^ixuiAo)^-!]^ for e = 
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Then the matrix Il^vC'^) admits the representation 

[A'/2] h [NjJl-h I \m-r\p I c\\\^^ 

For x-large enough, the matrices Y^^'^'^\X) andA]\f(A) depend smoothly on x and holomorphically on Ae U. One 

has the decomposition 

(C.3) nj'-*'^^ U) - J] [e (V,; e)] ^ nj!;^'^^ {X) [e (v,; e)]"? where v^^±q and vq = ^ . 

ee{±l,0) 

r/ie matrix Y)^^'^'^^ {A) is such that it does not contain any oscillating factor in its entries. Moreover, for all k 
there exists an N-independent constant C > such that 

(C.4) \\An\\l«(u) <— and 

These estimates also hold for the first order partial derivatives ( in respect to x or A). 



n 



{m,b,p) 
N:€ 



< C^x*^^ with w = 2 max ] sup |!R (y - v (e^))| 



Proof — 

Recall that the matrix 11// can be represented in terms of Cauchy transforms (or their + boundary values) on 
(C.5) n^v iA)^C^^o...o C^^ [I2] (A) = U) . 

Above and in the following, [M] (A) for A ^ corresponds to the case where in the integral representation 
(6.3) for this operator we substitute the + boundary value with A ^ 'rf. This is clearly a well defined expression. 
We decompose the jump contour for IT according to £n - dD U £n with dD - dDq^s U dD-q^s U dD^^g. 

The exponentially small in x terms gathered in A^f can be written as 

(C.6) An (A) = [Csof''''' ° {c~] ° {C^S^h] U) , 

whereas 

(C.7) (A) - An (A) = V C^^ ° • • • ° C^^, . 

There, the sum runs through e e S^f = {e = (ei, . . . , e^) '■ e {±1,0}}. One can readily convince oneself that 
for any matrix function M such that AM e (Zuj there exists a constant c' such that 

(C.8) C~ [M]\\ <c'||AM||,2^ ^ . 

Thence setting c = max |c', c(Zn) , c(5£))| (we recall that for a curve F, c(r) stands for the norm of the + 
boundary value of the Cauchy operator on (F)), one gets 

\\An\\lo^(U) < ^^^J^^^^^^sn^^) With N^(A) = IIAIIi^^) + IIAIIi..^^) . 
Thus, the claim follows for A^r as, by construction, N^^ (A) = O {x'°°). 



K.K. KOZLOWSKI 



57 



It remains to obtain estimates for the remaining, algebraically small in x, part. For this we set 



°"3 



A (z) = [e (if (z) ; e (z))]"- A (z) [e (<p (z) ; e (z))] " with 



<p(z) - 
eiz) = 



Then, by carrying out similar expansions to (B.14) and (B.16) it is easy to convince oneself that 

2 

(C.9) ' ''"''"'^ ^ 



Above Aflf, stands for the ab entry of A. Also, there appears A instead of A as the oscillating factors have been 
already pulled-out, as in (C.2)-(C.3). Also the primes ' in front of the two sums are there to indicate that these are 
constrained. Namely, one should sum-up only over thoses choices of e e £/v and ka ^ {1,2}, a = 2, ... ,N which, 
upon the replacement A i-^ A would give rise to the oscillating factor associated with [nl?'^'^^], . By using 

the continuity of the -i- boundary value Cauchy operator on (dD), one shows that there exists a constant c such 
that for any e, T e {±1, 0} and any / e (dD): 



(C.IO) 
Then, 

(C.ll) 



n 



(m,b,p) 
N;e 



N-l 



< i^£)!:L7v,2>(A) . 



Since there exists c' > such that || A \\L'^(aD)< c'x^, the claim follows. Also, we stress that, by construction, 
njy!^^'^^ does not contain any oscillating terms in x in its asymptotic expansion when x — > -l-oo. ■ 
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